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The modified Hansen theory for the motion of an artificial satellite, as developed by Peter Musen, 
represents the gravitational perturbations. The development takes the form of several trigonometric series 
with numerical coefficients and specified arguments. The result of the modified Hansen theory is the position 
vector of the satellite which is obtained by evaluating these series. 

In the present paper the velocity vector and the six classical osculating elements are shown to be obtainable 


also from the series in the general theory. 


The ability to obtain osculating elements and/or velocity from the modified Hansen theory is considered 
a significant adjunct to that theory. Comparison with other general perturbations theories and with special 
perturbations techniques is now possible. Accurate transition from the general theory to numerical integra- 
tion makes it possible to transfer to the latter method when additional disturbing forces, not included in the 


general theory, become of consequence. 


INTRODUCTION 


HE modified Hansen theory (Musen 1959) 
represents the gravitational perturbations of an 
ificial satellite. The development takes the form of 
everal trigonometric series with numerical coefficients 
nd specified arguments. The result of the modified 
Jansen theory is the position vector of the satellite 
yhich is obtained by evaluating these series. The 
resent paper contains the derivation of the velocity 
rector and the osculating elements from this theory. 
comparison of the results of the modified Hansen 
heory with those of other general perturbations 
eories and with those of numerical integration is 


ow possible. 


_ DETERMINATION OF THE VECTORIAL CONSTANTS 
FOR THE EQUATOR 


A basic trihedron, (P,Q,R), connected with the 
sculating plane is represented by the vector toward 
igee, P, one directed perpendicular to the orbit, R, 
ind a third, Q, at right angles to P and R, and forming 
ith them a right-handed system. The components of 
ese unit vectors are, in terms of the Keplerian 


elements: 
+cosi sinw cos6-+ cosw sin (1) 
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Taking the following relations into consideration, 
(Musen 1959) which represent a modification of 
Hansen’s basic relations (Hansen 1862): 
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w=xX-C 
a=a0— (a—n) AE—N—K 
= 0)— (at+n)AE—N+K, 

and using the “mean” or Hansen values of w and @, 
(w)=m0—o0t (yta—n)AE (4) 
(8)=00— (a+n)AE, (5) 


the developments for the osculating values of P, Q, 
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and R are easily obtained. As an example of the result 
we give the development for P-: 


P,,=-+cos*4i cos2K cos[x—m0— yAE+ (w)+ (6) J 
—cos*3i sin2K sin[x—mo— yAE+ (w)+ (6) J 
+ sin?4i cos2N cos[x—ao— yAE+ (w)— (6) J 
— sin’47 sin2N sin[x—m—yAE+ (w)— (6) J. 


From Eq. (16) of the ““Application of Hansen’s Theory,” 
(Musen 1959) the osculating values of the \ parameters 
are 

\i=sing7 cosV d3=cossi sinK 


\o=sin3z sinN 4=cossi cosK, 


and their numerical values are obtained by evaluating 
their respective series. Defining 


B=ecos(x—mo—yAE) (6) 


y=e sin(x—m—yABZ), (7) 


we see that 


eP.=B{ (AP—A3*) cos[(w)+ (6) ]—2Asdz sin[ (w)+ (4) ] 
+ (APv—)2?) cos[(w)— (8) J— 22 sin[ (w)— (6) J} 


—v{(A2—)3") sin[ (w)-+ (6) J+ 2A3A4 cos[(w)+ (6) ] 
+ (Av—)¢2?) sin[_(w) = (0) J+2r1d5 cosL_(w) — (6) ]}. 


We also obtain similar expressions for the y and zg 
components of osculating eP and for the components 
of eQ. We then define (P) and (Q) to be the coefficients 
of 8 in the expressions for eP and eQ, respectively. This 
is done for two reasons. (1) We know x to be the angle 
from X, the departure point, to the osculating perigee 
and mo+yAE to be the same angle for the auxiliary 
ellipse to the perigee of the auxiliary ellipse. Hence 
x—mo—yAE is small and of the order of the perturba- 
tions, @ is of the order of e, y is approximately zero and 
(P) approaches the real value of P. (2) We note the 
following relationships exist: 


eP=8(P)+7(Q) (8) 
eQ=6(Q)—7(P). (9) 
We define (R)=R since R does not contain @ or y. 


One of the results of the Hansen theory is a 
trigonometric series, 


W=X+Y sinF+Z cosF, 


where X comprises those terms in W independent of F, 
and Y and Z are composed of the terms in W which have 
a coefficient for F of +1 in the argument. We can 
obtain Y and Z easily by evaluating the series W—X 
at F=0 for Y, and at F=}r for Z. Using the expressions 
for Y and Z (Hansen 1838) 
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we can easily obtain 8 and y. 
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Finally, in order to obtain the osculating values of 
Q and R, we have 


(P:)=+Ae—Az) cosl(w)+ (@) ] 
— 2dsd4 sin[(w)+ (@)] 

+ (Ar—Az*) cosl(w) — (6) J 
— 21d» sin[ (w)— (4) | 
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: + 2 (AzA\a+A1A3) cos (w) 

(Qz)= — (A2—A3) sin[(w)+ (6) ] 
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+212 sin[ (w)— (6) ] 
(Qz) = +2(Ai\s— ee) cos (w)— 2(Agds-F Aids) antos 
(Rx) = +2(Aa+Asd3z) sin (@) —2(Ads—A1A3) cos (6) 
(R,) =) (AqAa+Asd3) Cos (@) —2 (AgAa— Ada) sin (@) y 
(R.) —) Ag+A3— —)7). ( 
It is interesting to note that the above nine values 
(P), (Q) and (R) are nothing more than the n 
components of the I'-matrix (Musen 1959) which 


used to rotate the auxiliary fF from the auxiliary ( 
(Q), (R) system to the inertial system. 


(P2) (Qz) (Re) 
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DETERMINATION OF VELOCITY 


Le 
I 


For the determination of the osculating values 
the velocity vector from Hansen’s theory we start v 
the familiar Laplacian integral 


RX v+hr+heP=0, ( 


where v is the velocity vector, r° is a unit vector in 
direction of r, and h=[a(1—e?) }-’. Since the posit 
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ctor of the auxiliary satellite f is in the same direction 
leetse., 
r(t)= (1+»)F(z), 


ere (1+) is merely a lengthening or shortening 
ctor, then 


=P = (P) cosf+ (Q) sinf. 


RXv+h(P) cosf-+h(Q) sinf-+heP=0. (17) 


7 a vector cross multiplication of R=(R) into the 
yove equation we obtain osculating v in terms of 
), (Q), and Q; namely, 


v=h(Q) cosf—h(P) sinf-+heQ. 


ais equation represents the form the Hamiltonian 
ntegral” takes in the Hansen theory. From Eqs. (9), 
|2), and (13) we find 


1 
heQ= heo(Q) + ——{ (1—e')(Q)¥— (P)Z}. (19) 
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re velocity then becomes 
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0 
‘yuo @ now from Eqs. (15) and (20), 
i sinE hk 
wicl a 
lary 1—e cosE hyo 
“y= (22) 


T: (1—e7)} cosE h oe : 


[ao]? 
1—e) cosE ho 2 


THE OSCULATING ELEMENTS 


The determination of the osculating elements is an 
‘tra result since one could obtain them from r and y. 
Owever, it is felt that the theoretical expressions for 
ie elements in terms of the constants of integration, 


, €o, and zp and some of the resulting series of the 
\eory, Aa, Az, As, As, W, ho/h, would be useful. We 
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- Finally from the osculating value of the eccentric 


anomaly, obtained in the equation above, we find 
M = Eose—e Sin E oso. 


We have applied the formulae derived in this paper 
to the Vanguard I satellite for a particular moment of 
time. Using a computer program coded by Mrs. Jaylee 
Burley, we have found that the osculating elements 
obtained from the position and velocity vectors are in 
complete agreement with those obtained from the 
immediately foregoing formulae. 


CONCLUSION 


The ability to obtain osculating elements and/or 
velocity from the modified Hansen theory represents a 
significant adjunct to that theory. Comparison with 
other general perturbations theories and with special 
perturbations techniques is now possible. Accurate 
transition from the general theory to numerical inte- 
gration makes it possible to transfer to the latter 
method when additional disturbing forces, not included 
in the general theory, become of consequence. 
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In the numerical integration for orbits, it has been proposed that the secular variation of the orbit 
elements be computed in multirevolution steps. The present paper develops an algorithm for a k-revolution 
step which is analogous to the algorithm for a step of integration by Adams’ method. The analogue of the ( 
derivative generated at each integration step by Adams’ method is the increment to the orbit elements for f 


one revolution generated at each k-revolution step. The coefficients in Adams’ integration formula are 


replaced by polynomials in 1/k. These polynomials are the same as -those-appearing as coefficients in vi 
Lubbock’s summation formula and the constants in these polynomials are the Adams coefficients. 


1. INTRODUCTION 


N the numerical integration of differential equations 

governing orbital motion it has been proposed 
(Taratynova 1957; Thomas 1958) that the change in 
the orbit elements for one revolution be determined by 
conventional numerical integration over one revolution 
and that this secular change in turn serve as the basis 
for a multirevolution integration step. Taratynova 
treats the secular change for one revolution as a 
derivative with respect to revolution number and then 
numerically integrates with respect to revolution 
number. Our method, while similar to that of 
Taratynova, treats the secular change for one revolution 
as a finite difference, and introduces a new Adams-type 
corrector formula which reduces to Adams’ formula 
only in a limiting case. The same refinement to the 
Taratynova method has been independently developed 
by Mace and Thomas (1960) and numerically verified 
by them. In their study, they applied a multirevolution 
predictor formula to the first and second sums of the 
acceleration components, rather than a corrector 
formula to the orbit elements. Our work further 
supplements theirs through the development of a 


general formula [Eq. (15) ] for the coefficients in the — 


formula for the multirevolution increment, which we 
evaluate and present through the tenth difference 
terms. 

Our multirevolution step is:started by extrapolating 
the orbit elements & revolutions ahead. Then, starting 
with these extrapolated values, the differential equa- 
tions of motion are numerically integrated over one 
revolution. This integration yields the increment to 
the orbit elements for one revolution at the end of the 
k-revolution step. The extrapolated values of the orbit 
elements are then improved by use of our corrector 
formula, which is in terms of this increment for one 
revolution and its backward differences with respect 
to the k-revolution spacing. As in Adams’ method, the 
corrector formula may also be applied to an extra- 
polated difference table to obtain the extrapolated 
orbit elements required at the beginning of the step. 
In the next section the corrector formula is derived. 


In the final section the growth of roundoff erro. 
discussed. 


2. DERIVATION OF CORRECTOR FORMULA 


Let y, be an orbit element evaluated at the enc 
the 2th revolution; let the increment for one revolui 
be 


fn=Ynti—Yny 


and let the backward differences with respect to 
k-revolution spacing be defined by 


Vin=fn—fnoky n=k, 2k, 3k,: Pt. 
Further let factorial polynomials be defined by 


xen) = 


TT (“=F r). 


By the Newton-Gregory formula for interpola’ 
using backward differences, 


Summing on i yields the k-revolution increment in 
orbit element 


—yo= dD fri 
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We now wish to develop further the formula for 
coefficients A. One has 


Ao=1 
A= —3L1+ (1/k)] 
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re factorial polynomial in Eq. (8) may be expanded 


0 a power series 
iXV@ 5 ive 

() -Es(-), ( 
k r=l k 

vere S,’ is the Stirling number of the first kind 


irdan 1939a). These numbers satisfy the recurrence 
ation 


Sq Ss," J Si (10) 
rting from 

| So=1, (11) 
1 

| SYP=S7=0 for r>0. (12) 


Next, in Eq. (8) after substituting from Eq. (9), 
| sum on 7 may be expanded into a power series in k 
Bernoulli’s formula (Jordan 1939b). 


r+1)\ Bog 
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‘re the sum on q is to the greatest integer in 7/2 and 
| Bernoulli numbers By, are those beginning 


(Bo, Bo, By --)=(A, 1/6, =/30,-* 2 
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Jpon expanding Eq. (8) into a power series, first in 
| by Eq. (9) and then in 1/k by Eq. (13), one finds 


(—)? are i OS? sr+1 
im re, B De ( )\- eel (19) 
gq! 0 m2ar-+1\ 2¢ ka 


te the 1/k term to be expected from Eq. (13) does 
appear because it contains 


i 

a; 0, 7>1, (16) 
r=1 
1 factor. 
Che X’s through Ayo are thus found, by Eqs. (6), (7) 
| (15), to be the following polynomials in 1/2. 
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3. RELATIONS TO OTHER FORMULAE 


Equation (5) is related to Lubbock’s summation 
formula, the coefficients \; being essentially the same 
in both formulae (Steffensen 1927). The Lubbock 
formula, however, is more complicated and moreover 
would not serve our purpose. 

The interval defining f in Eq. (5) is an orbit revolution 
only in our particular application. Actually Eq. (5) is 
a valid solution for more general systems of difference 
equations, and the interval /# defining f may be 
arbitrarily chosen. As h approaches zero with hk held 
fixed, the limit of Eq. (5) is the increment to y expressed 
in terms of differences of its derivative. The limit is, 
therefore, the Adams corrector formula. One finds, 
upon multiplying the right member of Eq. (5) by h/h 


f 
y(hk)—y(0)= lim > hk); ( ve) 
h0  j=0 z=hk 


h 
hk fixed (28) 
dy 
SIH, gimye(= ~) : 
i=0 k> Gs mh 
Thus ; 
(=): j Sj 
lim \j= DD (29) 
Roe j! mor+l 


are the Adams coefficients. From inspection of Eqs. 
(17)—(27) the limit as k approaches infinity may indeed 
be recognized as the Adams coefficients. 

The foregoing is an analytic presentation of multistep 
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Adams’ Integrotion 


4 zy 


Multistep Summetion (nck) h nh 


Fic. 1. Comparison of Adams’ integration with 
multistep summation. 


summation and its relation to Adams’ integration. The 


corresponding geometry is shown in Fig. 1. In the 
Adams part of the figure the curve is a polynomial 
fitted to dy/dx at intervals kh, and the shaded area is 


nh a) 
— } dx=y,—Vn-t- 
fess: (- 


In the multistep part of the figure, the curve is a poly- 
nomial fitted to f/h at intervals kh (k=3) and the 
shaded area is 


(E)nin 


If kh is held fixed but & approaches infinity, f/k ap- 
proaches dy/dx and it is apparent from the figure also 
that multistep summation reduces to Adams integration. 


(30) 


(31) 
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4. ACCUMULATION OF ROUNDOFF ERROR 


In the multi-interval stepping, the growth of ror 
off errors can be adverse. In one multi-interval 


one has 
Vy, 
and, taking variances, 
| ovo. 
On the other hand, by & single-interval steps, 


k-1 
Vo fz: 
= 


Assuming that each f; has roughly the same vari 
and neglecting the correlation of the roundoff e 


in successive f’s, 
oy ko. , 


which is less than the variance in the multi-interval s 
for fixed variance of f. The variance of f can, howe 
be considerably reduced by summing the increment 
the orbit elements during the numerical integra 
over one revolution rather than by differencing 
orbit elements as implied by Eq. (1). The roun 
error of f thus becomes a small fraction of the varia: 
of y over a revolution rather than the same fraction 
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the possibility that dust and gas may exist in a significant 
ount in globular clusters is considered. It is shown that the 
y method of detecting such intraglobular matter at present is 
its obscuring effects, especially if the material has collected 
) clouds within the cluster. Evidence is presented that such 
aids exist in globular clusters, appearing as dark regions and 
es in the cluster. It is found that statistical fluctuations in the 
jected stellar distribution of the cluster cannot account for 
se dark regions. If these obscuring effects are attributed to 
ground clouds we must expect at least twenty such clouds 
square degree of sky at high galactic latitudes. It is concluded 
t such an explanation cannot account for all the obscured 
jions that have been found in globular clusters. 


I. INTRODUCTION 


“LOBULAR clusters are thought to be free of 
J interstellar dust and gas. It is important to 
umine this conclusion critically because of the 
nificant position ascribed to globular clusters in our 
trent evolutionary concepts. Such a consideration 
ms the subject matter of the present paper. It will 
shown that interstellar matter may indeed exist in 
ibular clusters; the observational evidence leading 
this conclusion is presented in Sec. II. A possible 
irce of intraglobular matter and processes which 
id to remove such matter from clusters are discussed 
Sec. III. The effects of intraglobular material on 
» stellar population of a cluster together with 
idence that such effects may exist are presented in 
2. IV. A summary is given in the last section. 


| Il. OBSERVATIONAL DATA 


Interstellar matter may be detected by two general 
ithods: by its obscuring effects and, under the proper 
jiditions, by its emission characteristics. The appli- 
‘ion of these methods to globular clusters is considered 
the present section. 


(a) Obscuring Matter 


Evidence that dark lanes and regions exist in 
‘bular clusters is at least 100 years old, dating back 
the Earle of Rosse’s (1861) visual observations of 
13. He published a remarkable sketch showing a 
‘ies of dark lanes within the cluster. Interest in this 
yblem was renewed about the turn of the century 
Holden (1891) who, using photographic methods, 
© called attention to the dark lanes in M13. An 
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The source of the intraglobular clouds can be satisfactorily 
accounted for by the mass shed by stars evolving to the white- 
dwarf stage. About 0.5% of the total mass of the cluster becomes 
available in this manner in the time between passages of the 
cluster through the galactic plane. 

The effects of such intraglobular matter on the stellar population 
of the cluster are discussed. Such clouds may condense into 
stars and thus account for those stars which lie in the vicinity of 
the main sequence and above the turnoff point in the color- 
magnitude diagrams of some globular clusters. These stars appear 
to have a nonrandom distribution in the cluster, supporting the 
above suggestion as to their origin. 


interesting summary of these and other similar 
observations is given by Miss Clerke (1903). 

More recent observations of obscured regions in 
globular clusters are those of FitzGerald (1955) who 
comments on some dark lanes in w Centauri, and 
those of Idlis and Nikolsky (1959) who discuss M4. 
In the latter case the authors consider the possibility 
of statistical fluctuations in the projected stellar 
distribution of the cluster as a possible explanation of 
the dark regions in M4. They find that it is highly 
unlikely that such fluctuations would account for the 
size and number of such lanes and conclude that there . 
is obscuring matter in the cluster. 

Mrs. Hogg (1959a,b) has also drawn attention to 
the dark regions in globular clusters. She, however, 
suggests that these obscuring effects may be due to 
(presumably close) foreground clouds, pointing out 
that globular clusters present an ideal background 
against which such clouds would be detected. 

To examine this problem in more detail all globular 
clusters with galactic latitude >15° and lying within 
the declination range of the Palomar Sky Survey 
(5>—27°) were examined. Two sources were used, the 
Sky Survey prints and the Lick-Crossley plate file. 
These sources complement each other; the Crossley 
plates with their greater scale and brighter limiting 
magnitude allow an examination of the more central 
region of the cluster whereas the Survey prints, reaching 
a fainter limiting magnitude, yield information on the 
outer region of the cluster. 

Although factors such as image growth, exposure, 
and cluster concentration class affect general appearance 
of a cluster, it is possible to identify regions in a 
number of clusters which have a high likelihood of 
being caused by obscuration. Thirty-two globular 
clusters lie within the specified declination and galactic 
latitude limits. At least 12 of these clusters show one 
or more prominent dark regions and lanes on either or 
both the Crossley plates and Sky Survey prints. These 
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TaBLeE I. Globular clusters containing dark regions and lanes. 


NGC M NGC M NGC M 
1904 79 5904 5 6341 92 
2419 6171 6402 14 
5024 53 6205 13 6864 75 
$272 3 6254 10 7078 15 


clusters are listed in Table I and several are illustrated 
in Figs. 1 and 2. 

The 32 clusters examined in this survey cover 0.8 
square degree of sky. This figure is based on Mowbray’s 
(1946) estimates of their diameters. If all the regions 
designated as obscured are to be attributed to fore- 
ground clouds and not to intraglobular matter or 
statistical fluctuations in the distribution of stars 
within the cluster, then the number of such foreground 
clouds per square degree of sky, and at the higher 
galactic latitudes considered here, will be 1.25, where 
N is the number of obscured regions identified in the 
32 clusters studied. A lower limit for V is 16. We 
would therefore expect at least 20 interstellar clouds 
per square degree of sky. Because their appearance is 
so pronounced, these clouds must be of rather high 
opacity. Neglecting for the moment the question of 
statistical fluctuations, we are forced to accept either 
of two conclusions: (1) Relatively opaque interstellar 
clouds at high galactic latitudes are much more common 
than heretofore thought, or (2) globular clusters 
contain interstellar clouds. 

On the basis of their general appearance and opacity, 
we may compare the foreground clouds proposed in 
(1) to globules. An extensive study of the frequency 
and distribution of globules has not been made. The 
preliminary survey of selected areas in the region of 
the Milky Way by Bok and Reilly (1947) will, however, 
serve as an adequate basis for comparison. By searching 
in areas of background nebulosity or high star density 
they concluded that globules ‘‘seem to occur at the 
rate of one per square degree.” At least twenty times 
this number is needed to account for the obscured 
regions found in globular clusters. This large difference 
makes (1), above, unlikely. Further, such a high 
projected density of foreground clouds would predict 
obscuration in at least some elliptical nebulae. This is 
not the case. Galaxies such as NGC 205 and NGC 5128 
(radio source 13S4A) which are classified as peculiar 
cannot owe their peculiar appearance to foreground 
clouds; the unusual features are obviously associated 
with the galaxy itself. 

The obscured regions are prominent in photographs 
taken over a wide range of exposure time. A significant 
feature in the longer exposed photographs is the 
frequent appearance of stars superposed on an obscured 
region. It should be noted that even with the appear- 
ance of these stars in obscured regions there often 
remains a significantly large area where there is a 
striking absence of stars. From star counts in the 
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general field, it is immediately evident that these sta 
cannot be field stars in front of a noncluster obscurin 
cloud; the number of stars per square degree of sky” 
far too low. It has been suggested (Hogg 1959a) that th 
appearance of these stars represent cluster stars shinin 
through a nearby foreground cloud, such stars becomin 
more visible with longer exposure. For such a situatio 
one ‘would require a large fluctuation of opacity in th 
nearby small obscuring clouds. In addition there wo 
have to be a cluster star, of sufficient brightness tok 
visible, located in exactly the same line of sight as f 
low-opacity region of the nearby foreground clout 
Such a set of circumstances seems very unlikely. | 
would be more reasonable to expect that longe 
exposures would bring up faint cluster stars which al 
situated in front of an obscured region located withi 
the cluster. 

Thus, if we are to attribute the obscured regions t 
foreground clouds we require: 


(i) The occurrence of such clouds to be twent 
times greater than that of globules in the galactic plam 
(ii) An explanation of their absence as seen again: 
other bright backgrounds such as elliptical nebulae. 
(iii) Large fluctuations in the opacity of such cloud 
so that we may see cluster stars through the clouds. ' 


Or, we may attribute the obscured regions to intré 
globular clouds. Since the necessary mass for suc 
clouds can éasily be accounted for, in fact is expecte 
(Sec. III), we find it reasonable to conclude thé 
interstellar matter, concentrated into clouds, exists 1 
globular clusters. { 

(b) Statistical Fluctuations 5 

In the preceding discussion it was assumed that th 
obscured regions were not caused by statistic 
fluctuations in the projected stellar distribution of 
cluster. We will now proceed to show that such 4 
assumption is reasonable, i.e., that, on the basis ¢ 
statistical fluctuations in the stellar distribution, th 
probability of finding as large a star-free area, and E 
close to the cluster center as the obscured regions, 
vanishingly small. 

The most direct approach to such a problem is 
construct a two-dimensional model of a globular clusté 
containing the various parameters which describe 
cluster and see if relatively large areas are free ¢ 
stars. In theory an infinite or at least a very larg 
number of such models would have to be constructe 
before some definite statement as to the probabilit 
of finding large star-free regions could be mad 
However, there is another approach which may kt 
used. In this case only one cluster model (or even pai 
of a model) need be constructed. From this model w 
may obtain the frequency distribution of differer 
sized star-free regions. Such a distribution may b 
expected to approximate a normal distribution (fc 
positive values only). From the sample data we ma 


fic. 1. Examples of globular clusters showing obscured regions near the outer part of the cluster. NGC 6171 was photographed 


‘h the 48-inch Schmidt (Copyright National Geographic Society—Palomar Observatory Sky Survey). The M13 photographs were 
cen with the Crossley Reflector, the blue by N. U. Mayall, the red by A. G. Mowbray. Note particularly the regions at position 


gles 5° and 210° in NGC 6171 and at 130° in M13. 


Fic. 2. Examples of globular clusters showing obscured regions near the outer part of the cluster. Both clusters were photograph 
with the 48-inch Schmidt (Copyright National Geographic Society—Palomar Observatory Sky Survey). Note particularly the regi 
at position angle 0° in M10 and 215° in M14. 
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ien evaluate the dispersion of this distribution and 
lus assign a probability to the finding of any sized 
ar-free area. It is this second method which was 
nployed here. 

The following information is necessary for the 
mstruction of a globular cluster model: 


(1) The projected stellar density gradient, N(R). 

(2) The cluster luminosity function, ¢(M). 

(3) The limiting magnitude, mm, of the model. 

(4) The relation between image size, d, and apparent 
agnitude, m. 


model patterned after M13 was constructed using 
le following input data. Shapley’s (1915) star counts 
_ the cluster gave N(R,m)ARAm, the number of 
ars brighter than magnitude m in Am and within 
mcentric annuli of size R4AR. A cluster luminosity 
imction (M3, Sandage 1957) allowed extrapolation 
om Shapley’s magnitude limit to m~18.5 (photo- 
aphic). The Greenwich formula 


m=a-+-bd? 


as used to relate image diameter, d, to apparent 
agnitude; this approximate relation is adequate for 
\€ present purpose. The constants a and 0b were 
‘rived from measures on Mt. Wilson 60’’ photographs 
‘M13. From these data a table of V(REAR, m+Am) 
its equivalent V(R+AR, d+Ad) was constructed. 
random number drawn from Triplett’s table (1952) 
id representing the polar coordinates of an individual 
ar was then assigned to each of the WN stars. The 
sulting model for two quadrants of the annulus 
SR<2’ is shown in Fig. 3. 

‘The frequency of different sized star-free regions in 
is model is shown by the histogram in Fig. 4. The 
ita are sufficiently well represented by a normal 
Stribution with a standard deviation o=61; the fit is 
own in the figure. This particular annular region for 
'13 was chosen because it contains a very prominent 
itk lane visible at position angle 130° in Fig. 1. The 
€a of the lane in a photograph of the above limiting 
agnitude is about 460 square seconds of arc; it is 
(own as an X in the frequency distribution of Fig. 4. 
a the basis of this model, this obscured region is 
tger than 70; the probability of finding such a large 
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Fic. 3. Model globular cluster. One-half of the annulus 
1’<R<2’' is shown. 
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Fic. 4. Frequency distribution of different sized star-free 
regions in the model globular cluster. A normal curve with c=61 
is fitted to the histogram. The X represents the size of the obscured 
region in M13. 


lane strictly on the basis of statistical fluctuation is 
essentially zero. Several of the obscured regions in the 
clusters listed in Table I are comparable in size to the 
M13 lane and a similar conclusion may be drawn as 
to their reality. 

It should be pointed out that a cluster model con- 
structed in the preceding fashion overestimates the sizes 
of star-free regions. This results from a failure in the 
model to account for the turbidity effect in the photo- 
graphic emulsion. For the highly crowded stars in a 
globular cluster this overestimation may be quite large. 
This lends further support to the conclusion that the 
more centrally located regions in globular clusters 
which are being considered here are not due to statistical 
fluctuations of the stellar distribution. 

In summary we may say that the dark lanes and 
regions in globular clusters are very likely real; i.e., 
they are not due to statistical fluctuations. They may 
be due to the obscuring effects of nearby clouds, but 
on the basis of their frequency and of the appearance 
of “foreground” stars it is more likely that these 
clouds are actually located within the cluster. The 
lack of a mechanism by which such clouds might form 
cannot be used as an argument against their existence, 
for then we face a similar problem as regards the 
interstellar clouds in the general field. 


(c) Mass of Intraglobular Clouds 


We may assign a lower limit to the mass of an 
intraglobular cloud from knowledge of its dimensions 
and opacity by assuming that the observed obscuring 
effects of the cloud are produced wholly by dust 
particles of optimum size for scattering and absorption. 
Adopting the values employed by Bok (1948) in 
assigning masses to globules we obtain 


App =6.7X 10°56, (1) 
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where A», is the absorption for photographic light in 
magnitudes per kiloparsec and 6 is the dust-particle 
density in g/cm’. We will use as input data values 
derived for the dark nebula in M13; these values are 
representative of the obscured regions in a number of 
globular clusters. The total photographic absorption is 
at least five magnitudes; a conservative estimate of 
the size of the projected image of the cloud is 
~1pcxX3 pe. We will assume that the cloud is a 
cylinder seen perpendicular to its major axis. The 
cloud is more likely to be spherical, the observed 
elongated form being caused by images of neighboring 
and foreground bright stars. We will thus under- 
estimate the total mass of the cloud by assuming that 
it is a cylinder of the above dimensions. 

The resultant dust-particle density is 6=7.5x10~* 
g/cm’ and the total mass of the dust particles 
Maust/IMo > 2.6. Adopting Wy1/Mause= 10, we obtain 
Metoua/IMo > 260. There may be several clouds in a 
cluster; only those lying in that part of the cluster 
closest to the sun would be detectable; the total mass 
of interstellar matter within a globular cluster would 
then be proportionally higher. A value of ~10*o 
would be a reasonable estimate of the total mass of 
interstellar matter in a globular cluster. 


(d) Emission from Possible H1 and Hi Regions 
in Globular Clusters 


We have assumed in the above calculation that 
most of the intraglobular matter is hydrogen. If this 
hydrogen is neutral and atomic it should be detectable 
by its 21-cm radiation. A search for such radiation 
from the globular clusters M3 and M13 was made with 
the 85-foot radio telescope of the National Radio 
Astronomy Observatory (Roberts 1959). A conserva- 
tive upper limit (i.e., it may be well below‘ this limit) 
to the antenna temperature for both of these clusters 
is 0.4°K. For small optical depths, 7, this antenna 
temperature yields upper limits to the mass of H1; for 
M3 it is 7009 and for M13, 200Mo. 

These results must be viewed with caution. Although 
they are based on an observational upper limit to the 
antenna temperature, they are also dependent upon the 
assumed opacity. 

The assumption of small 7 for such clouds may be 
incorrect; this together with an expected low value of 
the A kinetic temperature in a globular cluster would 
result in an underestimation of the amount of Hr 
present. Further, a significant amount of hydrogen in 
a molecular form may be present in the obscured 
regions of globular clusters and would thus be un- 
detectable. The formation of molecular hydrogen 
would be favored in the highly opaque and presumably 
cool regions considered here. 

Our current lack of knowledge of the opacity, gas 
kinetic temperature, and the ratio of atomic to molecu- 
lar hydrogen make an interpretation of the above 
21-cm antenna temperature uncertain. 
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Globular clusters are known to contain stars whis 
have relatively high surface temperatures; these 4 
the stars on the blue part of the color-magnitu 
horizontal branch. These stars may be expected - 
form significant Hit regions when surrounded - 
interstellar matter. A characteristic emission line fre 
such Hi regions is 43727. This line has not be 
observed in the integrated spectra of globular cluste 
(Mayall 1946). It therefore is important to con 
the expected intensity of this line to see if its absen 
is consistent with the concept of intraglobular matt 

Minkowski and Osterbrock (1959) have consider 
a similar problem, the occurrence of 3727 in t 
spectra of some elliptical nebulae. Lacking data + 
individual exciting stars in these nebulae they assum 
that the relative frequency of such stars is the same 
in globular clusters. They then computed the expect 
intensity of \3727 for the globular cluster M3, for whi 
the necessary data were available and applied tl 
result to their discussion of E nebulae. Their rest 
will be adopted here, but with an important modifi 
tion which will be noted below. 

We wish to compute the total expected emission 
\3727 in ergs/sec and compare this value to the ener 
contained in a corresponding wavelength region in t 
continuous spectrum of a globular cluster. Minkows 
and Osterbrock have evaluated the emission ¢ 
efficient j,3727 in ergs/sec cm’, for an Hm region wi 
an electron temperature of 10000°K and an electr 
density less than 10*/cm*; they obtained j)3727=1. 
X 10-7? ergs/sec cm’, where Ny is the densit 
hydrogen ions or atoms per cubic centimeter. Th 
have also computed the ionizing radius for the 
cluster stars which would make a significant cont 
tion to the ionized volume which, with the clust 
luminosity function, yields the total ionized volun 
This together with jx3727 gives the total energy emitt 
per second in \3727. 

Implicit in these calculations is the assumption th 
the interstellar matter is uniformly distributed throm 
out the cluster. Only then would the full ioni: 
effect of each star contribute to the total A3727 ae 
We have seen, however, that intraglobular matter r 
be primarily concentrated in clouds. For such 
situation the ionizing effect of an exciting star situa 
at some distance from the cloud will be reduced by 
factor proportional to the square of the distan 
between star and cloud (Strémgren 1948). Special 
if S is the ionizing radius for the case of a star loca 
in the cloud, then for a star at a distance 7 from t 
cloud (r large compared to the dimensions of the clou 
the extent of ionization within the cloud R is given | 


R=S8/39. ( 


Minkowski and Osterbrock tabulate individ 
values for S*Vq? for each star which will make 
significant contribution to the total ionized volu 
The sum of these values (corrected for incompleaal 


1 the cluster area studied) is }S°Vq?=4.5X10°. We 
fill adopt this value of 2S*NVy? as a measure of the 
mizing effect due to one star at a distance r from one 
loud. In general this will overestimate the emission 
rom Hi regions. The product, E3727 of the volume of 
ne ionized region and j)3727 is 8.1 10*/r® Vz! ergs/sec; 
ne ionized volume is assumed to be a sphere of radius 
'/2. Eys727 is evaluated for various values of r and Nz 
1 Table II. Only those entries for which r>cloud 
idius are tabulated. 

These values of E3727 are to be compared to the 
xpected continuous emission from the cluster in this 
favelength region. For this we will again make use 
f the calculations of Minkowski and Osterbrock. They 
nd that the energy emitted by M3 near 5540 is 
'ss40= 1.710% ergs/sec A. From the relative energy 
istribution of globular clusters (Roberts 1956) as 
erived from the six-color observations of Stebbins and 
Vhitford (Stebbins 1950), we find E)3727/Eys540= 0.63 ; 
ence the energy emitted in the continuous spectrum 
ear 43727 is 1.110% ergs/sec A. Comparing this 
alue with those of Table II we find that only if the 
loud ion density is very low or the exciting star very 
lose to the cloud would we expect to detect \3727 in 
mission. From the very high opacity of the cloud we 
tould expect Vy>>10. We conclude that \3727 would 
ot be observable in the spectra of globular clusters. 
. similar calculation for Ha emission from Ht regions 
1 globular clusters yields the same conclusion. 


Ill, FORMATION AND REMOVAL OF 
INTRAGLOBULAR MATTER 


The source of intraglobular matter must be the 
luster itself. The high velocity of globular clusters 
‘lative to interstellar clouds makes it virtually 
npossible for a cluster as a whole to pick up inter- 
‘ellar matter by an accretion process. However, there 
‘Treason to expect a significant source of intraglobular 
latter from those cluster stars which are evolving to 
xe white-dwarf stage and which, in the process, shed 
art of their mass. These considerations have pre- 
lously been discussed by Takebe and Matsunami 
1957) and by Sandage (1957). The latter author has 
stimated the total amount of matter which has been 
jected into intraglobular space in the lifetime of the 
uster. He obtains a value of ~10°o for M3, which 
about one-half the present mass of the cluster. This 
laterial is not expected to remain in the cluster. 
assages by the cluster through the galactic plane will 
e very effective in removing the intraglobular matter. 


TaBLe II. Emission of \3727 in ergs/sec for various 
values of r and Vz. 


\Nu/cms 
aS 10 102 5X10? 
3 1.7102 
} 5 5.0 10% 8.210% 
10 8.0X 1055 8.010% 1.310% 
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If globular clusters follow highly eccentric orbits 
about the galactic center (von Hoerner 1955), then 
M3 will spend about 4X108 years away from the 
galactic plane; most of this time the cluster will be 
close to apogalacticon. The average amount of inter- 
stellar matter in M3 is then 


1 10°X4X105 
~X—————-=2X10°Mo, (3) 
2 10” 


where the age of the cluster is taken as 10 years. 

The value found in Eq. (3) represents a time average 
over the life of the cluster. The amount of interstellar 
matter to be expected in a globular cluster for a recent 
epoch is given by 


ANX (AM/At) X37, (4) 


where AN is the total number of cluster stars which 
lie in a particular part of the color-magnitude diagram, 
ASM is the mass loss suffered by these stars during their 
evolution, and Aé is the time necessary for a star to 
evolve through the particular region being considered. 
T is the time between passages through the galactic 
plane. For that part of the color-magnitude diagram 
which lies above the main-sequence turnoff point 
AN«At or AN/At is constant. We may therefore 
evaluate Af for any convenient region of the C-M 
diagram, thus 


At=6.2X108XON i 


m1) 


222) dg, 
=, (5) 
L 
where X is the mass fraction of hydrogen in the star, L 
its luminosity and q; is the fraction of the total stellar 
mass located in the core (Hoyle and Schwarzschild 
1955). The integral may be numerically evaluated 
from the model sequence computed by Hoyle and 
Schwarzschild (1955, Table 8) for type II stars. Such 
an integration yields At~6X108 years for the range 
M,~+2.4 to M,~—3.0. From the luminosity function 
for M3 (Sandage 1957) we estimate AV, the number of 
cluster stars which lie in this magnitude range and 
along the giant branch, as 3X 10°. 

We now require only Adv. For the initial mass we 
adopt 1.29%. The final mass will be that of a white 
dwarf. It does not matter if the final mass over the 
range considered here is actually that of a white dwarf, 
since for some Af this mass loss will be realized and we 
then make use of the relationship AV/At=constant, 
and may thus compute it for any At. For a white-dwarf 
mass we adopt 0.69%. This is based on Greenstein’s 
(1958) statistical determination of 0.569% for white 
dwarfs and on direct mass determinations which are 
available for three white dwarfs in visual binaries; 
these yield 1.0, 0.4, and 0.49M%o. For }T of 2X 108 years 
we would then expect 0.6X 10*Mo of interstellar matter 
in M3. Within the accuracy of the input data we may 
estimate the mass of interstellar matter in M3 as 
~10°Mo. 
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TABLE IIIa. The mass, I0/MMo, of a cloud which would contract 
under its own gravitational attraction for various values of 
Ny and T°K, p=1.4. 


ak 
Na/cm? 10 25 50 100 
10 500 2000 5600 16 000 
102 160 630 1800 5000 
5X 10? 70 280 800 2200 


Because of the low escape velocity from a globular 
cluster, ~10 km/sec, evaporation of interstellar gas 
from the edge of the cluster would occur at relatively 
low kinetic temperatures (S50°K). If such evaporation 
were to go on at a rapid rate, the above values for the 
interstellar mass in M3 would be far too large. However, 
if such material were to gather into clouds the evapo- 
ration rate would be greatly reduced. It was shown in 


Sec. II that there is good evidence that such clouds: 


(of dust, at least, and presumably gas) do exist in 
clusters and we therefore assume that evaporation is 
not the dominant mechanism for the removal of intra- 
globular matter. 

In summary, we may expect to find ~10°Mo of 
interstellar matter in M3. This is approximately 0.5% 
of the present mass of the cluster. Within the range of 
cluster masses and varying values of the semi-major 
axes of their orbits, this percentage may be assumed 
to be representative of globular clusters in general. 


IV. EFFECTS ON THE STELLAR POPULATION 
OF A CLUSTER 


Globular clusters have been described as a “dying” 
population—one in which stars evolve to their final 
state as white dwarfs or escape completely from the 
cluster. The presence of interstellar: clouds in the 
cluster may alter this picture. If these clouds are 
dense and massive enough, stars may form from this 
interstellar matter and thus add a small number of 
new stars to the total cluster population. 

We may employ the virial theorem to make a rough 
estimate of the minimum mass of a cloud which will 
contract under its own gravitational attraction. For 
contraction we require the second time derivative of 
the moment of inertia to be negative. This condition 
will be met if 

2E~<—i, 


where £ is the total kinetic energy of the particles and 
Q is the total potential energy. We shall assume that 


TABLE IIIb. Radii in parsecs of a spherical cloud of the 


above mass and density, .=1.4. 
“BK 
Nu/cmé 


10 29 50 100 

10 7.1 113 15.9 2268 
10? DES) 3.5 5.0 fii 
5X 10? 0.5 1.6 Die Sat 
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the initial distribution of material is roughly spheri 
and of radius R, temperature T, and density p. \ 
find that for contraction 


nII4 TNA TANS 
—>1.2x10-»(=) (-) : { 
Io u p : p 


where yu is the mean molecular weight. For a hydrog 
density Ny and y=1.4, Eq. (6) becomes j 


i 


I7/aMlo > 50(T!/Ny?). f 


Table IIIa contains values of 9%/Me for varic 
combinations of Ny and T. 

From previous considerations we expect the mass 
the cloud to be IT<10*M%e. We see from Table 0 
that we require relatively high densities and | 
temperatures. The former are not too extren 
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Fic. 5. Composite color-magnitude diagram for the blue st 
in nine globular clusters. The filled circles represent the 
stars, the open circles, the stars in the remaining eight clust 
The age-zero main sequence and the region of the M3 dwarfs ; 
subgiants is also shown. 


especially in view of the high opacity noted — 
intraglobular clouds. A low temperature may also 
expected to exist in such opaque clouds. 

Table IIIb contains the corresponding radii, 
parsecs, for a spherical configuration of mass a 
density corresponding to those of Table IIIa. We ag: 
find that, to be consistent with the apparent sizes 
the dark nebulae found in globular clusters, we requ 
densities of several hundreds and temperatures 
several tens of degrees. A cloud possessing a cent 
condensation rather than the above assumed unifo 
density would have a larger gravitational self-ener; 
The requirements on the temperature and dens 
would then be less extreme (Kahn 1960). 

The foregoing considerations apply to the entire clo 
We may then expect stars to form within this cloud 
a manner similar to that operating within clot 
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isting in the general field. In essence, we would 
pect a cluster to form within the boundary of a 
ybular cluster. If such a star-formation process had 
curred relatively recently we would expect to find 


(1) Main-sequence stars which are young compared 
the age of the cluster, i.e., stars which sit above the 
ister main-sequence turnoff point. 

(2) A clustering of such stars within the projected 
2a of the globular cluster. 


lere is evidence that both of these occur in globular 
isters. 

The occurrence of stars above the main-sequence 
moff point in the color-magnitude diagram of globular 
isters has been noted by Sandage (1953) and by Arp 
955). These stars, which are distinct from the 
rizontal branch stars, are too.common to be 
ributed to field stars superposed on the cluster. An 
cellent example of such stars is to be found in the 
lor-magnitude diagram of M3 (Sandage 1953, 
hnson and Sandage 1956). For convenience in the 
lowing discussion these stars will be referred to as 
ie stars. A composite color-magnitude diagram for 
e blue stars in nine globular clusters is shown in 
g. 5. The filled circles represent M3 blue stars while 
e open circles represent the blue stars in the remain- 
zy eight globular clusters. The source material for 
is figure is given in Table IV. All observations not 
eady on the B, V system were transformed to this 
stem except those for M53 which are on the P, V 
stem and close enough to the B, V system (Cuffey 
58) for the present purpose. The M53 star plotted 
B—V (really P—V)=—0.7 is an exception to the 
sumed similarity of the two systems. The age-zero 
4in sequence and the M3 main-sequence and sub- 
unt regions are also shown in this figure. The ap- 
rent separation into two groups of blue stars implied 


TaBLeE IV. References for data used in Fig. 5. 


No. blue Apparent B—V 
aster = stars modulus reddening Reference 
M3 31 15.67 0.00 Sandage, A. R. 1953, 
Astron. J. 58, 61. 
110 4 14.82 0.35 Arps Hay ©. 11955; 
Astron. J. 60, 317. 
[2 2 16.15 0.04 Arp,) Ey (Ge) 1955; 
Astron. J. 60, 317. 
[22 1 14.00 0.44 Arp, H. C. and Mel- 
bourne, W. G. 1959, 
Astron. J. 64, 28. 
[GC 3 16.85 0.10 Sandage, A. R. and 
1147 “Walker, M. 1955, 
Astron. J. 60, 230. 
5 2 15.08 0.14 Arp, H. C. 1955, 
Astron. J. 60, 317. 
115 1 15.89 0.00 Ary lie \G@.2/ 1955; 
Astron. J. 60, 317 
113 3 14.67 0.14 Acpsia Ha|Ca, 1955, 
Astron. J. 60, 317. 
[53 7 16.90 0.00 Gutteyw 9.2) 1958, 
Astrophys. J. 128, 
219. 
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Fic. 6. (a) Solid histogram: the observed angular distribution 
of the blue stars in M3. Dashed histogram: the blue star distribu- 
tion corrected for different sized samples of stars available for 
each sector. (b) Angular distribution of all stars with mp ,<19.00 
in Sandage’s catalogue for M3. These data were used to obtain 
the dashed histogram in (a). 


by the color-magnitude diagram is very likely due to a 
selection effect; many of the clusters were studied 
only to intrinsically bright magnitudes. 

It has been suggested (Sandage 1953) that such blue 
stars represent some small fraction of the total cluster 
population which, in their evolutionary processes, have 
remained well mixed and are therefore evolving 
essentially up the main sequence rather than to the 
right into the giant area. Such mixing would, presum- 
ably, be due to a very high rotational velocity of the 
star thus establishing the necessary meridianal circu- 
lation currents. However, such high rotational velocities 
are very rare for the type of star considered here and 
we must consider this suggestion with some reservation. 

A number of globular clusters were examined for the 
possibility of a clustering tendency among the blue 
stars. The results for M3 are shown in Fig. 6. This 
cluster has a good sample of such blue stars which are 
distinctly separated from the scatter of the main- 
sequence and subgiant parts of the color-magnitude 
diagram as well as from the blue horizontal branch 
stars. There are 29 such blue stars that (i) lie in a 
region above the turnoff point defined by the following 
limits : mp»<19.00, C.I.<0.38 and to the left of a line 
located by the coordinates m,,= 18.20, C.I.=0.50 and 
Mpy= 19.00, C.I.=0.26, and (ii) are listed in Sandage’s 
(1953) catalogue under Series I and II (Series III does 
not cover the entire cluster and therefore could not 
be used). 
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The distribution of these 29 blue stars within six 
equally spaced sectors is shown by the solid histogram 
in Fig. 6(a). Eleven of these blue stars lie in one sector. 
The number of all stars brighter than m,,=19.00 in 
Sandage’s catalogue is not evenly distributed among 
the six sectors. Except for crowding and variables, the 
stars included in the catalogue were randomly chosen 
(Sandage 1959) and do not necessarily indicate the 
actual angular distribution of the cluster stars. An 
incompleteness factor must thus be applied to the 
numbers of blue stars per sector to adjust for this 
uneven catalogue coverage. When normalized to the 
average number of all stars per sector, the number of 
blue stars becomes 30.5; thus, one would expect 
5.1+(5.1)? blue stars per sector. Instead, almost one- 
half of this normalized number of blue stars (13) are 
found in just one sector. The formal probability of 
finding 13 or more stars in a sector is less than 0.014. 
The normalized distribution of the blue stars is shown 
by the dashed histogram in Fig. 6(a). The distribution 
of all catalogue stars brighter than m,,=19.00 is 
shown in Fig. 6(b). 

Unfortunately, the other globular clusters for which 
C-M diagrams and catalogues are available show too 
few blue stars for any significant conclusion to be 
drawn. One possible exception is M10. This cluster has 
been studied by Arp (1955) to a magnitude limit only 
a bit fainter than the horizontal branch. He found 
four rather bright blue stars. These are numbers I-32, 
J-33, I-34 and III-48, and may be identified in Arp’s 
Fig. 3. Three of these stars fall within a circle of 18” 
radius; this is an area only 0.25% of the total area 
studied by Arp. If they are at the same distance within 
the cluster, they are separated by less than one parsec. 
These numbers must be viewed with great caution. 
Not only are the stars few, but other ,clusters which 
contain a similar small number of known blue stars 
do not show this clustering effect. Further, these 
bright blue stars may have an origin different from the 
fainter blue stars found nearer the main-sequence 
turnoff point. 

If these blue stars are indeed newly formed within 
the globular cluster we would expect that they would 
initially move as a group through the cluster. It is 
likely that such motion would carry the group close 
to the central regions of the cluster. The effect of this 
would be twofold: (1) the shorter relaxation time in 
this region would cause the group of stars to disperse 
and (2) any remnants of the cloud from which the 
stars were formed would also be dispersed. This latter 
mechanism can be expected to be operative for any 
intraglobular cloud in general. Since several oscillations 
through the cluster would be necessary for significant 
tidal disruption, the time scale would be similar to 
that for the period of oscillation through the galactic 
plane and so need not affect the previous arguments. 

The formation of the blue stars has to satisfy another 
requirement—that of the time necessary for the stars 
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to contract from some low-density condensation with 
the cloud. We may set as an upper limit for this tir 
the period of oscillation through a globular cluster 
a star (or cloud) situated near the edge of the cluste 
This is about 3107 years. The contraction time f 
the latest type star that might be considered here, F 
is ~1.3X10" years. Thus, a blue star may form in 
time short compared to the period of oscillation of 
cloud through the cluster. The presence and appare 
clustering of the-blue stars are thus an expect 
by-product of the interstellar matter in globular cluste! 

Another mechanism for creating and clustering t 
blue stars is accretion. The possibility that accretis 
could explain the occurrence of early type stars in t 
general field has generally been abandoned. T 
problem of the relative motion of stars and clou 
proved too formidable. In a globular cluster, tk 
problem does not exist. The velocity dispersion mu 


be relatively small, otherwise the star or cloud wou 


escape. It is thus important to examine accreti 
quantitatively to determine if it is a competing mech 
nism with star formation. 
The accretion rate (Bondi and Hoyle 1944) is giv 
by 
dM /dt =a (2rpG2M?/u'), ( 


where 9%=mass of star 
p=density of cloud 
u=relative velocity of star and cloud 
G= gravitational constant, 


and a@ is a numerical factor which depends upon t 
way in which the accretion stream has been establishe 
It has a value between 1 and 2. 

The cloud will retard the motion of the star; f 
u=0, Bondi (1952) obtains in place of Eq. (8) 

dM /dt=a' (2rpG?9n?/C*), ( 
where C is the velocity of sound at ea al 
0.5<a’ <2.24. 

The time, 4, and distance 2, for m, the initi 


velocity, to be reduced to a value somewhat less, th 
C have been derived by McCrea (1933). He obtains 


Uo? (1 
(4a+38) 2rpG?3o 
Uo" 
ak en (1 
(Sa+48) 2apG?IMo 


here, 9% is the initial mass of the accreting ste 
Following McCrea we shall adopt a=2 and B= 10. 

For a star trapped in the cloud (w<C) integratic 

of Eq. (9) gives the time for the star to increase 1 

mass by an amount ASM to some new value 1, thus 
CG et 

a= |e (1 
a! 2arpG? SC OM 


tg— 
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is then the total time for the process, i.e., to stop the 
ir and to add mass ASW. 
The conditions of the problem require that 


(i) «; be small compared to the cluster diameter, 
(ii) t2 be less than the time between passages through 
2 galactic plane, 

(iii) fg be less than the time for the blue stars to 
olve off the main sequence. 


r (i) we may take x to be of the order of one-tenth 
a cluster diameter, or about five parsecs. Values of 
to~10* years will satisfy (ii) and in most cases (iii). 
Tables Va and Vb give values of x; and /, for various 
drogen densities, Nw, and percentage mass increase. 
y evaluate Eqs. (10)—(12) the following were adopted: 
=1 km/sec, a’=2.24 (isothermal case), MWo=1Mo, pu 
1e mean molecular weight) = 1.4, and C=0.55 km/sec. 
The sequence of events which would have to occur 
uld involve the capture of a number of stars of 
out solar mass by a very dense and relatively large 
ud. Initially, these stars would have just a small 
locity difference from that of the cloud—enough to 
ing them into the cloud. Once in the cloud these 
irs would accrete varying and relatively small 
1ounts of mass and as a result increase their lumi- 
sity and temperature, thus moving them into another 
rt of the color-magnitude diagram. Eventually the 
ud would be destroyed or removed from the cluster 
s in the case of the formation of the stars) leaving a 
all group of stars whose kinematic and physical 
operties have been altered. As this group moves 
rough the cluster, interactions with the general field 
cluster stars will cause the group to disperse (again 
in the formation case). Thus both the occurrence 
d clustering of the blue stars can be accounted for. 
The data of Tables Va and Vb indicate that, within 
@ approximations of the input data, accretion can 
<e place. However, the requirements here (Vy and 21) 
2 more extreme than those for the case of regular 
w formation and we would expect that if we had such 
uds they would just as likely condense into stars. 
e may therefore conclude that accretion, though 
ssible, is not the principal mechanism for the forma- 
m of the blue stars. 

Throughout this discussion we have assumed that 
e admittedly scanty data on the clustering of the 
le stars represent a real clustering and not just a 
ance occurrence. From the data presently available 


TABLE Va. The time f2 (in years) for various values of Ny 
and the final mass of the star, I%/Mo. 


ee 

Niz/cmé ita Eales} 2.0 5.0 10.0 
10 fpommelO os wetos2 820.2 ~22.2)><108 
10? iecmeO-sm ido.2, 20.2 22.2 107 

m »<10 1S) 2.0 3.0 4.1 4.4 X10? 
108 (eco sa etoe2, 2052) 2232 <X10° 
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TaBLE Vb. The distance 1; and the corresponding mass of 
a spherical cloud of this diameter. 


Ny/cm? xpc I/We 
10 4.8X 10? 1.9X107 
10? 4.8X10 1.9105 
5X10? 9.6 7.8X 108 
108 4.8 1.9 108 


such an assumption seems reasonable. An obvious 
desideratum is more information on the numbers and 
distribution of blue stars in clusters; this will primarily 
come from color-magnitude diagrams which extend to 
relatively faint magnitudes. We need not expect that 
all globular clusters will contain blue stars. It also 
follows that if the reality of the clustering of the blue 
stars is firmly established, the possible explanations of 
the existence of these blue stars in globular clusters 
seems pretty well restricted to the two given here. It 
should also be kept in mind that there may be at 
least two distinct classes of blue stars, one possibly 
being represented by the subluminous O-type star in 
M3 and the other by the large number of fainter blue 
stars which are also found in M3. 


V. SUMMARY 


Various methods of detecting interstellar matter in 
globular clusters have been discussed. It was shown 
that Hu regions and their characteristic emission lines 
are very unlikely in a globular cluster. Further, the 
interpretation of even weak 21-cm emission is difficult 
because of a lack of knowledge of the physical conditions 
of intraglobular matter, especially the ratio of atomic 
to molecular hydrogen. The only method available at 
present is that of the obscuring effect of interstellar 
dust. 

There is evidence that such dust (and presumably 
gas) clouds are present in globular clusters. This is 
shown by the large number of obscured regions found 
in globular clusters. It is difficult to attribute all of 
these obscured regions to relatively close foreground 
clouds; the number required by such an explanation is 
too large. In addition, such close and opaque clouds 
would need to have large fluctuations in opacity to 
account for the appearance on long-exposure photo- 
graphs of stars superposed on the obscured regions. 
The likelihood that these regions are due to a statistical 
fluctuation in the stellar distribution was shown to be 
small. 

This intraglobular material can be attributed to the 
shedding of mass by stars evolving to the white-dwarf 
stage. A significant amount of such material will collect 
in a globular cluster between passages of the cluster 
through the galactic plane. The collection of such 
matter into dense clouds could result in the formation 
of new stars. This would account for those stars which 
sit above the main-sequence turnoff point in the color- 
magnitude diagrams of some globular clusters. Such 
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an explanation would further predict that these stars 
would have a nonrandom distribution in the cluster. 
The scanty data thus far available indicate that this 
is indeed the case. 
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Measures on nonpolarized stars are reported for the McDonald 82-inch and 36-inch and for the Goethe 


Link 36-inch reflectors. For different telescopes the amount of instrumental polarization differs greatly. 
The polarization increases toward the ultraviolet. The instrumental effects may cause photometric as 
well as polarimetric errors, and a method to avoid them (Thiessen and Broglia 1959) is referred to. 
Measures on depolarization effects are reported and small corrections are derived for incomplete action 
of the depolarizer disks and of the Wollaston prism. Depolarization by the mirrors of the 82-inch telescope 


is found to be negligible. 


INTRODUCTION 


4 


HE main purpose of this series of papers is to 
compare interplanetary and interstellar particles, 
with measures of polarization and color over a wide 
range of wavelengths. In order to cover a wide range 
of wavelengths, and to depend as little as possible on 
weather conditions, a Wollaston photometer was built 
that uses a depolarizer for calibrations; the ordinary 
and extraordinary rays are measured with two multiplier 
phototubes by simultaneous integrations. The photom- 
eter has been described (Gehrels and Teska 1960); the 
effectiveness of the depolarizer was checked with each 
of the wide-band filters in use. 

Preliminary results have been published for stars, 
reflection nebulae, Venus, Mars, the moon, and the 
sunlit blue sky [diagrams and further references were 
presented (Gehrels 1960a) at a Polarization Conference 
sponsored by the National Science Foundation and by 
the Lowell Observatory ]. The detailed observations will 
be given in this series. Paper I describes the corrections 
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for spurious polarization, introduced by aluminiz 
mirrors, and a method to avoid the effects in the futu 

Dollfus (1957) observed the polarization by telescc 
mirrors, of normally incident light; the following ce 
clusions are drawn from his work. (1) Dollfus has r 
found any instrumental polarization on silvered mirro 
silvered either chemically or in vacuum. (2) Soi 
aluminized mirrors, aluminized in vacuum, show stro 
effects, while others do not. (3) The polarieanieg var 
very ees over the mirror surface. (4) The effe 
appear to be caused by the coating rather than by 1 
perfections of the mirror polish. 

‘Thiessen and Broglia (1959) are studying the effe 
extensively, and preliminary conclusions are dra’ 
from their experiments. (5) The instrumental pola 
zation is not caused by the aluminum evaporati 
being at an angle different from normal to the surfa 
For steep angles, however, the conclusion may be 4 
ferent (see later). (6) The instrumental polarization 
caused by the glow discharge, for final cleaning of t 
mirror, being at an angle different from normal to t 
surface. Presumably, the asymmetric potential or gli 
discharge establishes an electric field in a thin surfe 


INSTRUMENTAL POLARIZATION 


yer, such that the aluminum is aligned in a crystalline 
ucture. (7) The amount of polarization is independent 
the thickness of the coating; the orientation of polar- 
tion is with the electric vector maximum parallel 
the electrode. Presumably, the crystals are aligned 
rpendicular to the electrode and the polarization is 
differential absorption. (8) Not all mirrors with 
lique glow-discharge show instrumental polarization, 
d the observed amounts vary widely. The polarization 


parently occurs only if the crystalline dimensions: 


ve reached a critical value with respect to the wave- 
eth of light. 

When the angle of evaporation differs greatly from 
rmal to the surface, by more than 70°, the above con- 
sions (i) and (5) may no longer be correct. Reimer 
157) deposited metals very obliquely, in layers with 
ckness of about 0.04 micron, and he found strong 
larizations in the transmitted light. 

For polarization measures of greater precision than 
ut of the present program, effects of image formation 
the cathode and of instrumental circular polarization 
ve to be taken into account. These effects have been 
cussed by Behr (1960) and by Serkowski (1960), 
pectively. 


THE OBSERVATIONS 


Ihe instrumental polarization is observed on stars 
ut are chosen because of parallax greater than 0702 
d/or because no polarization was found by other 
servers. At least five stars per observing period are 
served for this purpose. A table of corrections is 
tained from the observations, after change of sign, 
the individual photometer orientations (180, 270, 
), 240, 120, and 210 deg), an example of which is 
ywn in Table II of the description of the photometer 
ehrels and Teska 1960). In order to avoid the sys- 
natic effect of instrumental polarization, all our 
servations are corrected with such tables. 

After least-squares solution of the values in the cor- 
tion tables, but with signs of the original observa- 
ns, the instrumental polarizations are listed in 
ble I. Wavelengths are given in terms of 1/), the 
erse of the effective wavelength in microns. The 

“effective wavelength” is used here for 


A= SS(A)M/SS(A)AA, 


ere S(\) is the response of the tubes to the part of 
ite light transmitted by the filter, according to 
nufacturers’ specifications. The percentage polari- 
ion is P=100(1—T»2)/(1+72), where J; and J» are 
intensity of electric vector maximum and minimum, 
pectively; and @ is the equatorial position angle a 

electric vector maximum. The results of the 
Donald 36-inch in Table I agree closely with the 
rections for instrumental polarization published by 
nr (1960). The measures at the Goethe Link Ob- 
vatory were made by R. E. Samuelson as a part of 
rogram on Venus. 
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TaBLE I. Observations of instrumental polarization. 


ete 1/r Period > PB: 6 JP 
Apr. 17-22 0.06 5°4 
30! 1.01 Apr. 23-30 0.12 135.3 % 
; July /Aug. 0.20 67.9 0.18 
[Apr. 1960 0.32 146.2 
82” 1.46 July/Aug. 0.13 81.0 OTs 
Oct. 1958 0.05 tee 
Apr. 17-22 0.19 26.5 
82’ 1.79 Apr. 23-30 0.06 Til S 0.09 
July /Aug. 0.07 168.5 
Apr. 1960 0.08 Sif 
82” 2.38 July /Aug. 0.12 164.4 0.12 
Apr. 17-22 0.48 oe 
” Apr. 23-30 0.44 IAS) 
o2 2.79 July Aue. 0.45 168.9 0.48 
Apr. 1960 0.56 178.2 
82” 3.08 July/Aug. 0.63 ae 0.63 
” May 0.12 WPS 
pg EOL December 0.17 141.2} ae 
36” 1.46 December 0.12 Sieh 0.12 
Dec. 1958 0.35: tee 
” February 0.36 see 
so 1.79 May 0.64: 12.8: 0.38 
December 0.27 18.1 
36” 2.38 December 0.98 26.2 0.98 
ip May Due 18.2:) 
“6 2.79 \ December 187 20.8% 1:96 
36”” 3.08 December 1.97 23.7 1.97 
Oct./Nov. 0.08 99.2 
GL 1.01 Dec. 1-16 0.34 106.0 0.18 
Feb. 1960 0.12 147.3 
Oct./Nov. 0.07 124.0 
GL 1.46 Dec, 1-16 0.05 99.1} 0.09 
Feb. 1960 0.15 118.1 
Oct./Nov. 0.16 109.3 
GL 1.79 +Dec. 1-16 0.22 110.9 0.14 
Feb. 1960 0.04 141.1 
Oct./Nov. 0.12 90.4 
GL 2.38 Dec. 1-16 0.10 96.4 0.13 
Feb. 1960 Oe 120.7 
Oct./Nov. 0.14 106.3 
GL 2.79 Dec. 1-16 0.15 89.5 0.11 
Feb. 1960 0.05 90.6 
Oct./Nov. 0.19 102.1 
GL 3.08 {Dec. 1-16 0.16 76.6 0.19 
Feb. 1960 0.21 105.6 
2 82/’=McDonald 82-inch, at Cassegrain; 36’ =McDonald 36-inch, at 


Cassegrain; and GL =Goethe Link 36-inch, at prime focus. 
b Observing periods in 1959, unless stated otherwise. 


The precision of the mean values in the correction 
tables is found from the repetition of the individual 
determinations, namely +0.0014 mag. (p.e.) 0.0004; 
after least-squares solution of the values at 6 orienta- 
tions the probable error in Table I is +0.0007 mag., or 
+0.03%+0.01. This is to be compared with the 
probable error found from the results on strongly 
polarized stars, in Paper II (Gehrels 1960b), namely 
+0.07%. The precision of the results in this program 
is little affected by instrumental polarization, especially 
when, in future work, nonpolarized stars are observed 
more often. 

The scatter of the individual results about their 
average, for each telescope and filter in Table I, gives 
a probable error of 0.04%+0.01. The agreement of the 
above numbers +0.03% and +0.04% indicates that 
the instrumental polarizations depend little, if at all, 


= Fic. 1. Percentage po- 

2 larization, introduced by 

aluminized mirrors, ver- 

sus the inverse of the 

| - effective wavelength in 

82 microns. For the Mc- 

Donald 36-inch and 82- 

inch and for the Goethe 
ie) Link telescope. 
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on age of the coating, on dust on the mirror, or on the 
alignment of mirror and photometer axes. The Goethe 
Link mirror had been re-aluminized immediately before 
the observing period of Oct./Nov., 1959. The 82-inch 
mirror had been washed prior to the run of April, 1960. 
For each observing session, the axis of the photometer 
is aligned again with that of the mirror(s). However, 
these conclusions are based mostly on the data for the 
82-inch and Goethe Link telescopes, whereas the large 
effects of the McDonald 36-inch may be more uncertain 
and different for stars of different colors (Behr 1960). 


WAVELENGTH DEPENDENCE 


Table I shows large differences between different 
telescopes, but the three reflectors have the low values 
at longer wavelengths in common. The large scatter of 
the position angles suggests that some of the small 
polarizations at longer wavelengths may not be real. 
The wavelength dependence for each telescope is shown 
in Fig. 1, with a line drawn through the averages of the 
last column in Table I. 

If the instrumental polarization is due to differential 
absorption by crystalline elements, such absorption 
would occur twice, namely upon incidence of the light 
and again after reflection by underlying parts of the 
coating (thickness usually about 0.1. micron). The 
crystal diameter, 2a, may perhaps be estimated with 
calculations of the extinction. by aligned metallic 
cylinders (van de Hulst 1957). The strongest polari- 
zation in van de Hulst’s Fig. 69 occurs at 2ra/A~0.3, 
so that for the mirror(s) of the McDonald 36-inch 
2a~0.02 microns, and possibly smaller for those of the 
82-inch. With observations at wavelengths shorter than 
3000 A, one might obtain the shape of Fig. 6 of Paper II 
(Gehrels 1960b). This explanation, however, is very 
tentative; multiple-reflection phenomena might explain 
the observed polarizations. Electron-microscopic studies 
of the mirror surfaces, as in Reimer’s work, with a 
determination of dimensions and alignment of the 
crystalline structure, would be decisive. 


THE AVOIDANCE OF INSTRUMENTAL POLARIZATION 


Instrumental polarization may cause small systematic 
errors in photometry as well as in polarimetry. The 
photometric errors, caused by the spurious polarization, 
in stellar work with the McDonald 36-inch are less 
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than 0.001 mag. But if a telescope with 5% insti 
mental polarization were to be used on an object wi 
20% polarization, a reflection nebula for instance, 
systematic error up to 0.01 mag. could appear. In 
dentally, a Newtonian telescope may have about 5 
instrumental polarization. It is especially important f 
high precision in polarimetry to avoid the instrument 
effects. y 

Telescope mirrors usually are aluminized with t 
evaporation normal, at least symmetrical, to t 
mirror surface. In addition, the glow discharge shot 
be normal to the surface. Even so, a small residt 
polarization may be found, the cause of which is bei 
studied at the Hamburg Observatory (Thiessen, pr 
sonal communication). Thiessen and Broglia are usi 
electrode disks that are, parallel to the surface, slight 
larger than the mirror; only in one case until now ¥ 
there a small amount of instrumental polarization le 
The Goethe Link mirror was aluminized with comple 
symmetry, by Leroy M. E. Clausing of Skokie, Illino 
and the results are good, as shown before. Mr. Clausi 
uses the glow discharge sparingly as he has observ 
(personal communication) that it causes a thin surf 
layer, presumably consisting of oxidized oil. This ¢ 
served surface layer may actually be the one para 
by Thiessen and Broglia. 


TESTS OF DEPOLARIZATION 


The effectiveness of the depolarizer is occasiona 
checked by rotating a piece of Polaroid HN 32 (Polaré 
HR for the infrared) in front of the photometer wh 
the depolarizer is in the beam. The rotation causes 
change of 1.2 units, of the nearly 100 units amplif 
deflection for one of the tubes; this means a spuric 
depolarization of 1.2 in what is supposedly 10 
polarization. The effect is probably due to nonperf 
orientation of the optical axes of the two depolari: 
disks. The checks are made with our combinations 
tubes and filters for 1/A=1.01, 1.46, 1.79, and 2. 
No wavelength dependence has been observed. 
1/A=2.38 and 3.08, the checks can be only qualitat® 
using the strong polarization of the sunlit blue sky a 
source, because the Polaroid is not effective at th 
wavelengths. The depolarizer breaks down when 
much narrower filter is used, such as half-stande 
thickness of Corning 2030 combined with a 1P21 tul 

The effectiveness of the depolarizer and the Wollast 
prism combined was checked with another Wollaston 
front of the photometer. The percentage polarization 
the ordinary and extraordinary rays was determin 
with the IBM 650 as usual, and the results are listed 
the first part of Table II. Allowing for the fact that t 
Wollaston prisms are involved, a depolarization eff 
is found of 2.5 per 97.5% observed polarization (3.6 ] 
96.4% observed polarization in the ultraviolet), 
which in each case 1.2% is due to the disks of the « 
polarizer. The amounts for the present photometer m 
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TaB1E II. Observed polarization of injected light. 


Source Location 1/r Jz 

: % 
Vollaston prism, ordinary Directly into ee ae 
1.79 96.5 
"ray photometer 279 94.2 
Vollaston prism, extra- Directly into ae me 
ordinary ray photometer 279 94.0 
01 77.9 
290°=35” 79 81.0 
79 86.3 
01 77.9 
295°-23"” 79 80.7 


Directly into 
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~I 
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‘ransmitted by stack of anes Ge 82.1 
nine glass plates 01 idk O 
DOr 79 80.9 

79 85.1 

01 77.6 

70°28" 79 77.4 

79 85.4 

Directly into rc i is 

photometer 2.79 86.0 


ctually be slightly smaller as the second Wollaston, 
hat was used for the beam, has two cracks. The de- 
larization effect is probably caused by multiple 
cattering inside the Wollaston prism; the wavelength 
lependence of the effect seems to confirm this ex- 
lanation. 

Depolarization by the two aluminized mirrors of the 
2-inch telescope was looked for with a 1-inch beam 
ransmitted by a stack of nine thin crown-glass plates, 
he angle of incidence being the Brewster angle (for 
ach color). In otherwise perfect darkness the beam was 
ent, parallel to the telescope axis, either incident on a 
mall portion of the primary mirror or directly into the 
hotometer. The measures were made in August 1960; 
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they are listed in the second part of Table II. The loca- 
tion of incidence on the 82-inch mirror is given for the 
telescope east of the pier and pointed toward the zenith. 
For instance, 290°-35’" means a spot at azimuth 
(counted from the north through the east) 290° and at 
a distance of 35 inches from the center of the mirror. 
The measures were actually made with the telescope in 
a more horizontal, more accessible, position. After 
reflection by the mirrors the polarizations are not less 
than when the beam goes directly into the photometer; 
the mirrors do not show any depolarization. The larger 
scatter of the observations that have the light directly 
into the photometer is due to the more difficult align- 
ment in that position. Incidentally, the observed 
amounts of polarization are appreciably larger than 
the predicted ones (see Jenkins and White 1957), 
61.6% at 1/A=1.01, 62.7% at 1/A=1.79, and 64.7 at 
1/A= 2.79. 

The conclusions on depolarization are summarized 
with the following example: When the observed amount 
of polarization is 10.00%, the correction for depolariza- 
tion should be +0.26 (+0.37 in the ultraviolet), of 
which 0.12 is due to the depolarizer and the remainder 
to scattered light within the Wollaston prism. Alumi- 
nized telescope mirrors, probably in general, do not 
cause an appreciable amount of depolarization. 
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The probable error of the polarization measures is 0.07%, or 0.0015 mag., found from observations 
made under various weather conditions. The interstellar polarization has a maximum at 6500 A, decreasing 
sharply toward longer and gradually toward shorter wavelengths. The wavelength dependence is quali- { 


tatively interpreted with computations (van de Hulst 1957) of the extinction by aligned cylinders. Particle 


sizes of 0.05-0.49 microns are derived for various refractive indices; the agreement of observations and (a 
calculations is equally good for metallic and dielectric particles. 


N the first part of this paper the performance of the 
photometer is checked by comparing observations 
on strongly polarized stars with results obtained by 
other observers. The second part deals with a remark- 
able wavelength dependence found for the interstellar 
polarization. This wavelength dependence had been 
qualitatively predicted by Davis and Greenstein in 1951 
and was discussed more recently by Davis (1959). The 
infrared part has been observed by Hiltner (Strémgren 
1956), while some decrease in the ultraviolet has been 
found by Behr (1959). 


THE OBSERVATIONS 


Eight stars were observed for the purpose of a com- 
parison with existing data. The observations are listed 
in Table I; they were made with the McDonald 82-inch, 
except for the ones of May 2, 1959, which were made 
with the McDonald 36-inch reflector. The stars are 
arranged in order of galactic longitude, the symbols 
have been explained in Paper I (Gehrels 1960), and 
p=P/46.05. Colons are used when the average of the 
residuals in the least-squares solution is greater than 
0.30%. : 

The precision is found from internal agreément of 
the repeated observations, giving half-weight to the 
values with a colon. The amount of polarization has a 
probable error of +0.07%, or 0.0015 mag. ; the probable 
error of the position angles is +2°7. Some of the 
measures were made through clouds and with poor 
seeing ; the results were not affected. The precision for a 
faint object is, of course, better when more light gets 
through. The over-all precision may therefore be slightly 
improved by observing under excellent conditions only. 

It is seen in the fifth column of Table I that some 
nights give position angles systematically different from 
those of other nights. The differences, although always 
less than 5°, seem to indicate a systematic error of 
unknown cause. 

The last three columns of Table I give the comparison 
with results of other observers, taken from a summary 
of references, effective wavelengths, and data, by 
Serkowski (1960). The agreement with the other ob- 
servers is well within the observational precisions. Cor- 
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rections for depolarization, as derived in Paper I, ha 
not yet been applied to the present results. 


WAVELENGTH DEPENDENCE 


The polarizations of Table I are repeated in Table 
after normalization to the average amounts of H 
183143 at 1/A=1.46 and 1.79. Also listed are the galact 
longitude and latitude with respect to the Lund po 
Until more stars have been observed, one may conclu 
that the values in the infrared are nearly the same. T. 
values in the ultraviolet may differ, but this, if re: 
is a very slight effect. Another, but incomplete a: 
therefore omitted, observation of HD 42379 casts dou 
on the high value of P for 1/A=2.79. There is obvio 
need for more observations, of stars at different distanc 
as well as at different longitudes, and with carel 
determinations of depolarization corrections and of t 
effective wavelength for each star and filter. 

The averages of the eight stars are in the botte 
line of Table II; values with a colon were given ha 
weight. The averages are plotted in Fig. 1; the cir 
radius indicates the probable error of each point, 
determined from the agreement of the eight stars. 
flat maximum of polarization occurs near \= 6500 
The line through 1/A=1.01 may actually be steep 
with the turnover at longer wavelength than that 
Fig. 1; narrow-band observations between 1/A=1. 
and 1.46 would be useful. 


INTERPRETATION 


The interpretation of Fig. 1 is based on calculatio 
of the extinction by long cylinders (perpendicular | 


Fic. 1. Observed 
centage polarization, : 
erage of eight stars, a 
function of the inve 
of the wavelength 
microns. 


OBSERVED 


yw S& a oo 30 
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TaBLE I. Observations, averages, and comparison with other observers. 


Single obs. 
Date UT Repeated obs. and averages Compared with 
) 1959 1/r 1p 6 P p 6 p 6 Observer 
Apr. 28 1.01 299 7192 é 
pr. ; ; . is 5 
July 31 1.01 2.78 64.8 2.89 07063 68 
July 30 1.46 fi i 4.04; 0,008 ~~ 65 
5 ae re ieee 4.08 0.089 66 O™110 62° Hall 
y 0.077 66 P. Smith 
July 30 2.38 a. ae 3.98 0.086 62 v. P. Smi 
Apr. 28 2.79 S101. 170-1 
July 31 2.79 2.98 65.8 3005 20-0607 68 
July 30 3.08 a er 3.07 0.067 63 
Apr. 28 1.01 O72 2 88.7 
July, 31 1.01 DSi 1 285-0 D162 = 00087) 89 
Aug. 6 1.01 B64. 92:5 
July 30 1.46 a ic 3.750 0.081. 86 
Apr. 28 1.79 Burt 19.0 
45 July 31 1.79 3.65 85.8 3.68 0.080 939 Baaee ogee y Ae 
Aug. 6 1.79 3.68 90-8 v. P. Smith 
July 30 2.38 a: one 32 2 e0.074 “> 86 
Apr. 28 2.79 2.67 5, 14 
July 31 2.79 2:80 86.1 Sie "0062 89 
Aug. 6 2.79 2.94 89.1 
July 30 3.08 Ma io 2.60 0.056 87 
Apr. 30 1.01 3 wi 4.21 0.091 0) 
Aug. 1 1.46 6.03 177.0 5.99 0.130 179 
Aug. 5 1.46 5.94 1.4 | 0.126 178 Behr 
Apr. 30 1.79 a ins 6.08 0.132 1 | 0.143 179. Hall 
43 Aug. 1 2.38 5.68 178.7 5.68 0.123 1 F130) 8179) Hilénee 
Aug. 5 2.38 5.67 3.1 0.147 0  ¥-P. Smith 
Apr. 30 2.79 oe Eee Sie 7 POsi12: 0 
Aug. 1 3.08 3.94: 178.3: 
Aug. 5 3.08 4.61: 22 @2e7 10.093. 0 
May 2 1.01 1.85 1.8 
July 31 1.01 2.011751 193 (0.042 178 
Aug. 1 1.46 2.64 179.5 
Aug. 5 1.46 2.93 21 25495 9200! 1 
May 2 1.79 2.72 179.0 2.73 0.059 178 0.062 2 Bebr 
July 31 1.79 273 177.1 erie ees 
78 Aug. 1 2.38 2.76. 179.1 4 all 
ae ee sae ee 2-72 0.059 0| 0.061 5 Hiltner 
May 2 2.79 2.31 s3 
July 31 2.79 2.40 178.0 DESO TOON. © v2 
Aug. 1 3.08 233 177.5 
Auge 5 3.08 2.19 8.2 gh A : 
Aug. 8 1.01 oe 7 241 70.046. 119°) 
Aug, 2 1.46 go 0.071) 114 0.074 115° °Behr 
Aug. 8 1.79 07 = 20-074) 1 115 0.072 115 Hall 
91 Aug. 2 2.38 A 2520.071 2 113 0.075 117 Hiltner 
Aug. 8 2.79 2.02 0.063. - 117 0.074 | 119 Serkowski 
Aug. 2 3.08 2.74: 0.060 ‘117: 
Apr. 28 1.01 1.92 0.042 169 peed eee ery 
JApr. 28 1.79 2.51 0.055 172 a 
84 Apr. 28 2.79 2.23 0.048 «172 0.062 170 Hiltner 
eee : : : fe 0.068 172 Serkowski 
Apr. 29 1.01 1.85 0.040 168 
79 Apr. 29 1.79 2:63 0.057 174 ee 
Apr. 29 2.79 2290: 2 OF065= 108: ‘ 
Apr. 29 1.01 2.00 0.043 179) 
17 Apr. 29 1.79 2:83 0.061 176 Hops ee emcee 
Apr. 29 2.79 2:23 0.048 178 a 


nce), made by van de Hulst and by Mrs. van P. 
th (van de Hulst 1957). In Figs. 2-6 the observed 
\rization is given with the drawn curve, copied from 
1 without change. The open circles give computed 
its as obtained from van de Hulst’s Figs. 69 and 71, 
circle radius indicates the precision of transfer from 
original figures. The difference (Qi—Qz») is taken, 
re (; and Q» are the efficiency factors for extinction 


when the electric vector is, respectively, parallel and 
perpendicular to the longest particle axis. (The com- 
parison of such points with the observed curve is 
rigorously valid only for a single size of particles, not 
for a distribution of sizes.) A factor f is chosen so that 
the maximum of /{(Q:—Q2) matches the maximum of 
the observed polarization. The circles in Figs. 2-6 give 
f(Qi—Qz2) plotted against 1/A=x/2mra, where x= 2ra/n 


472 


Fic. 2. The drawn line 
is copied from Fig. 1. 
The points with open 
circles are obtained from 
van de Hulst (1957), 
Fig. 71. 


Fic. 3. Same caption as 
for Fig. 2. 


Fic. 4. Same caption as 
for Fig. 2. 


is the variable used by van de Hulst, and a is the radius 
of the cylinder. The value of @ is chosen by trial and 
error to obtain the best fit. 

The results of curve fitting and the resulting values 
of 2a are in Table III. Van de Hulst (1955) fitted these 


TaBLeE II. Normalized percentage polarization and galactic 
coordinates. 


Normalized polarization at 1/d 


HD TRO UR RAG ee 22382 Ome on0S Be bt 
161056° ~4:19* 6515" 5292 S577 4.35) 454593462) =102 
154445 4.25 6.09 5.98 5.55 4.66 4.22 347 +21 
183143 4.21 5.99 6.08 5.68 5.18: 4.28: 21 —1 
198478 4.22 6.10 5.98 5.95 5.16 4.94 53 +1 

21291 3.90 6.02 6.04 6.01 5.40 5.06: 109 +4 
43384 4.62 --- 6.03 5 Ole 156 +5 
42379 4.24 6.03 6.65: 157 +3 
AGM, Ey ooo (i0R oni AL yG) doc 157 +1 
Average 4.24 6.07 6.01 5.79 5.09 4.57 ++» --- 
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Fic. 5. Same caption 
for Fig. 2. 
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Fic. 6. Same capti 
as for Fig. 2, but ir 
van de Hulst cS 


= 141-1 4li Fig. 69. 
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same extinction calculations to the observations 
interstellar reddening and he obtained the particle siz 
that are in the last column of Table III. The prese 
results appear consistent with those on interstel 
reddening. The cylinders are aligned, presumably | 
the galactic magnetic field. Van de Hulst (1950) h 
indicated that the results of the calculations remain t 
same, but for the amount of polarization, when the | 
cylinders are replaced by prolate spheroids. Mrs. vat 
Smith kindly sent me a few calculations made 
refractive index m=1.50—0.25 7, these also fit — 
observations and they give 2a=0.16. ; 
The fits appear equally close, within the uncerta: 
of the method; for different values of the refract 
index. It is therefore not possible from these data al 
to decide between dielectric particles, metallic 
ticles, or a mixture with suitable size distribution. 
mixture of particles will be considered in a followi 
paper, using the observations on reflection nebulae 


TaBLE III. Comparison of the observations with calculati 


(van de Hulst 1957) for cylinders. , E 
Figure 2a, from 2a, fr 
Refractive in this interstellar interste 
Substance index paper polarization reddeni 
Ices 125 7 0.49 
Impure ices 1.25-0.107 3 0.31 
Glass 1.50 4 0.28 
Impure glass 1.50-0.107 5 0.20 
1.50-0.257 tee 0.16 
Metallic 1.41-1.417 6 0.05 


NER Sek LAR 


| as the present results. However, the numerical 
cision of fitting the calculations to the peculiar 
ve of Fig. 1 is of the order of +0.02 in 2a. A choice 
particle sizes and refractive indices will be quite 
ited by the observed wavelength dependence of 
rstellar polarization. 

n conclusion, I would like to thank Mr. R. L. 
nfery for his help with reductions and illustrations. 


EF ASTRONOMICAL JOURNAL 


VOLUME 65, 


POLARIZATION 473 


REFERENCES 


Behr. A. 1959, Z. Astrophys. 47, 54. 

Davis, L. 1959, Z. Astrophys. 47, 59. 

Gehrels, T. 1960, Astron. J. 64, 466 (Paper I). 

Hulst, H. C. van de. 1950, Astrophys. J. 112, 1. 

. 1955, Mém. Soc. Roy. Sci. Liége, 4th Ser., 15, 403. 

——. 1957, Light Scattering by Small Particles (John Wiley & Sons, 
Inc., New York). 

Serkowski, K. 1960, Lowell Obs. Bull. 4, 318. 

Strémgren, B. 1956, Astron. J. 61, 45. 


NUMBER 8 OCTOBER, 1960 


Six-Color Photometry of Ten Classical Cepheids* 


Sotirios N. SVOLOPOULOS 
Lick Observatory, University of California, Mt. Hamilton, California 


(Received July 13, 1960) 


Six-color observations of the cepheids U Aql, FF Aql, SU Cas, X Cyg, SU Cyg, DT Cyg, Y Oph, S Sge, 
Y Sgr, and T Vul, made with the 22-inch Tauchmann reflector, are presented. Light and color variations were 
determined in each of the six colors. Values of the phase retardation of the various light curves are given. 
The relation between the light range and the inverse wavelength was found to be linear in all cases, but 
with deviating points in the colors U, V, and I; the deviations are attributed to absorption in the outer 
layers of the stars. Eggen’s type C cepheids are distinguished from those of types A and B of the same 
period in having a less steep slope in the relation between light range and inverse wavelength. Temperatures 
at maximum and minimum light have been derived from the six-color data: the average temperature range 
was found to be 1550° for type A and B cepheids and 750° for type C cepheids. Mean radii are given for 


four cepheids having accurate radial-velocity curves. 


URING 1957 and 1958 a program of six-color 
observations of ten bright classical cepheids was 
ried out with the Tauchmann 22-inch reflector of 
Lick Observatory. An infrared-sensitive Lallemand 
tomultiplier, designated L II by Kron (1958a) was 
d with a combination of filters that will be described 
1 future paper by Stebbins and Kron. Table I lists 
cepheids observed, in order of increasing period. 
2 elements were taken chiefly from the General 
alogue of Variable Stars (Kukarkin et al. 1958). The 
heid types given in the last column have been taken 
n Eggen, Gascoigne, and Burr (1957). 
‘he observations were made on dry, clear nights and 
standard extinction factors for each of the six colors 


TaBLE I. Cepheids observed in six colors. 


Epoch Period Vz Cepheid 
tar JD (days) (max.) (min.) AVz___ type 
‘Cas 2430404.134 1.949319 5.70 6.08 0.38 ? 
‘Cyg 24305.124 2.49934 5.60 5.90 0.30 C 
Cyg 33095.944 3.845664 6.39 7.10 0.71 A 
‘Vul 19372.151 4.435578 5.35 6:04 0.69 A 
Aql 24703.115 4.470959 5.13 5.48 0.35 C 
Sgr 2S. d 1335) -5.28° 6:08 0.80 A 
_Aqi 10170.173 7.023836 5.98 6.78 0.80 B 
' Sge 29091.48  8.3821723 5.24 5.99 0.75 B 
Cyg 21511.892 16.385680.. 5.82 6.89 1.07 B 
‘Oph 08694.86 17.11884 5.77 6.29 0.52 C 
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were used in the reductions. The observed colors of the 
cepheids were reduced to the standard six-color system 
by comparison with the six-color data for nonvariable 
stars given by Stebbins and Kron (1956). The com- 
parison stars used are listed in Table II. Two of these 
stars have been observed by Eggen (1955) in two colors, 
and their visual magnitudes in this system (Vz) were 
used in deriving visual magnitudes for the cepheids, by 
interpolation between green and violet. The observed 
visual magnitudes of the cepheids at maximum and 
minimum light are given in Table I. From 35 to 40 
observations in each of the six colors for each cepheid, 
a mean light curve was drawn for each color. From the 
differences between each individual observation and the 
mean of the nearly simultaneous B, G, and R observa- 
tions, a color curve was derived for each color. Figure 1 
shows these curves for SU Cygni; they may be taken 
as typical examples of the data derived in this inves- 
tigation. The mean deviations (in mag.) of the obser- 


TasLeE II. Comparison stars. 


Spectral 
Star HD type Vz 
15 Sagittae 190406 dGl 5.76 
35 Cygni 193370 F5 Ib 
HR 7807 194335 B2 Vp 
44 Cygni 195593 FS Iab 
¢ Cygni 197989 KO III 2.48 
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vations of SU Cygni from the mean curves.are listed in 
Table III. 

Table IV gives the brightness expressed in magni- 
tudes in each of the observed six colors, taken from the 
mean light curves, at decimal phases. The median 
magnitude between maximum and minimum was used 
as the zero point. Also given in the table are the 
maximum and minimum magnitude and their differ- 
ences. Table V gives the colors in the six-color system, 
as read from the smoothed color curves. The usual zero 
point, B+G+R=0, as previously stated, has been 
adopted. The slopes of the transformation lines for 
converting the observed colors to the standard system 
did not differ appreciably from unity. The extreme 
values were 1.07 for U and 0.91 for V. 

Inspection of the light curves reveals, usually, a 
phase difference (infrared being later) between U and I. 
In Table VI an estimate of this phase retardation is 
given for each of the ten cepheids. It is difficult to judge 
the accuracy of this quantity, since it is determined 
from light curves in which the time of maximum light 
is often uncertain owing to observational errors. From 


TaBLeE III. Mean deviations of SU Cygni observations (in mag.). 


U Vv B G R I 


+0.06 +0.03 +0.03 +0.02 +0.03 -+0.03 
=-0:05; »==0/02) 010% 722001 20:01 3:0:02 


Light curves: 
Color curves: 
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Fic, 1. Light and color curves 


for SU Cygni. 


the table, however, where the stars are listed in ord 
of increasing period, we see that there is no correlati 
between the phase difference and the period; also the 
is no correlation between the phase difference and t 
range of the mean light curve. ; 

Miss Hack (1956) has shown that from 4000 to 6000 
the light range of a cepheid is a linear function of t 
inverse wavelength in which it is measured. Plottt 
the ranges given in Table IV as a function of 
reciprocal of the effective wavelength, we find that t 
points for the three colors B, G, and R lie on a straig 
line for all the cepheids, while the U point deviat 
from this straight line for some of them, and the V al 
I points for all. Figure 2 shows a plot of this type f 
SU Cygni. : 

The deviation of U is negative, ie., the range 
magnitude is too small, presumably because of increas' 


Frc. 2. Light range 
a function of uw for 
Cygni. 
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TABLE IV. Six-color light curves of ten cepheids. 
hase U Vv B G R I Phase U V B G R I 
SU Cassiopeiae Y Sagittarii 
0P0 0.28 90.29 —0.21 —0.18 —0.16 —0.17 070 = 0.68 —0.62; 0144) —0'34_ —-0.26) —0,24 
0.1 = One On 1Sm-0n10- 0.14. —0.13. —0.17 Oot —0.62 —0.51 —0.38° —0.33 -—0.23 —0.18 
0.2 =-0,07))—0.05 —0.07 —0.05 —0.07 —0.13 0.2 M31 920; 19)" 0) 18:, 0.17, 9-0. 14 0.07 
0.3 +0.17 +0.08 +0.03 +0.03 +0.01 —0.08 0.3 —0.19 +0.03 0.00 —0.06 -0.07 —0.02 
0.4 20.25.9019 +0.12 +0:10 +0.07 —0;02 0.4 250505" 0.47 4-012 “-F0104" = 0500 » 0.04 
0.5 [0:2894-0.28 +0.20 -+0.16. +0.13 +0.04 0.5 40.27 40.33 +0.23 +0.14 +0.09 +0.11 
0.6 nog 0.29 022° 40:18 +0.14 +0.14 0.6 +0.48 +0.53 +0.37 40.24 +40.19 +0.18 
0.7 +0.24 +0.20 +0.15 +0.14 +0.09 +0.12 0.7 +0.73 +0.64 +0.46 +0.34 +0.27 +0.25 
0.8 40.24 40.06 +0.03 +0.06 +0.02 —0.02 0.8 +0.32 +0.43 +0.34 +0.27 +0.21 +0.19 
0.9 -0.0f 7-010 —0,f1 —0.08 —0.09 —0.11 0.9 —0.07 +0.03 +0.00 +0.01 +0.02 +0.06 
iL Mis us as ¥ at Max 72) 0564) 0,47 0536 2.0.27 0-005 
3 ee ae ee oe Min +0.73 +0.64 +0.46 +0.35 +0.27 +0.25 
Min 0.29//-+-0.30 0.22 +0.18 +0.14 +0.18 | ,,. ais 
Mix 0.59! 0.60 0.44 0.36 0.31 0.36 | Minax 28 0.93 0.71 0.54 0.50 
DT Cygni U Aquilae 
OPO —0,63, —0.63 —0.47 =0:39 0.28 = 
Eo nog 0s) 0:08 |g pas OSs Oso Loon cae wes 
0.1 gee 0200.16 0:13 0.11. 0.09 0.2 30: 203s 2017 2010 0s | L019 
0.2 =0,09' =0.12 —0.10 —0.08 —0.09 —0.06 : 5 Oo 1s e007) OO Oe ek 
0.3 E12) 0-11) 20,07+ 0:02 0108" —0-00 
eo +0.06 0.00 —0.02 —0.01 —0.04 —0.02 0.4 +0.03 +0.04 +0.03 0.00 —0.08 —0.08 
0.4 +0.19 +0.13 +0.08 +0.06 +0.02 +0.03 0.5 20.20 40.27 +0.20 40:18 +0:06 —0.01 
0.5 40.26 +0.24 40.16 +0.12 +0.09 +0.07 0.6 40.44 +0.52 +0.33 +0.27 +0.16 +0.09 
0.6 40.25 +0.23 +0.18 +0.15 +0.10 +0.09 0.7 Ee063 40.63. +0,45 +0:36 4-024 40.17 
0.7 +0.15 +0.10 +0.10 +0.10 +0.05 +0.06 0.8 BOL? 44-0149. 40:41 £0135, 1-027 40.22 
0.8 +0.02, —0.03 -—0.03 0.00 —0.01 +0.02 Ms 
i ee ein 0 00 20 07 2.0 08 0.9 40.08 —0.05 +0.01 +0.05 +0.07 +0.06 
; : . ; : Max O563 0/6340: 47 *— 0:39" 0) 2850008 
Max que OMmeeO 20m OS Oto, —Ol12; 0710 Min +0.63 +0.63 +0.47 +0.39 +0.28 +0.22 
Min. +0.26 +0.26 +0.18 +0.15 +0.11 +0.09 | Min-Max Te 20vne le 20reOFO4s NOs Sen cOb Oman Oss 
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TaBLeE V. Color variations of ten cepheids. 
Phase U V B c R I Phase U V B G 
SU Cassiopeiae Y Sagittarii 
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0.1 +0.37, —0.02 -+0.04 0.00 —0.04 =0)13 1) 0,140.50" 0,08" 40.07 n oo 
0.2 +0.41'°+0.03 40.05 0.00 —0.06 —0.47 | 012 40.52 -40/19) 2000, Seo 
0.3 +0.54 -+0.07 +0.07 +001 .—0.00 —020140:5° -+0.@" 40.28 015) oe 
0.4 +0.55 +0.11 +0:09 06.00 =0,10 —0.21 1 0.4 °40.73— --0:384-0,19)) om 
0.5  +40.55 £0.13 40,115 =000f =040 0.921 0s 408s 6 se 0.00 
0.6 +0.54 +0.15 40.144 -—0.01 =0:12 ©2007 1) 0:6" 0:08 6-48 --0.22) en 
0.7 +0.49 +0.13 40:12 (2001 —O.91- —0 201) oe) Gs 0.00 
0.8 +0.39 40.07 +0.07 0.00 ——0.07--=0.16 | 0/8 0.81" +0.40 +-0.1s) “Fos 
0.9 +0.36 —0.02 +0.04 0.00 —0.05-.—0:13| 0.9 40.61 40:22, 40:11) -tom 
DT Cygni U Aquilae 
Pe 
oFO 0.04 =0.28 =0.10 —0.01 40.12 +0.20] | 0*? Bie ae a ae 
0.1 =0.01 -0.25 -0.40, —0,01 (04) ROT 5: Gigge 4g a0) 0.20 0.00 ¢ 
0.2 +0.03 —0.22 —0.08 0.00 +0.08 +0.14] 9°; 14°93 10:47 40.28 0.00 ; ; 
013.) 0:00" 20:19 ~=0.07 0.00 +0.06 +0.11] 974 44712 40°53 10.31 6.00 : 1 
0.5 +0.15 -0.08 -0.04 0.00 +0.03 +0.07) 9°5 [1'53 40:71 40.34 0.00 =0.34 2mm 
0.6 +0.12 -0.10 -0.03 —0.01 +0.04 +0.06] 6°, 11°39 40:72. 40.34 0.00 2034) 2m 
0.7 +0.08 -0.14 -0.05 -0.02 +0.06 +0.12] 973 44°93 40.61 40.31 0°00 \igieaus aa 
0.8 +0.04 -0.21 —0.08 0.01 +0.10 +017! 9°) 14'or 40.41 40.25 0. 0OmmeeOnS “ai 
0.9. 0,00. 0.26, —0:90. | —Ov01 |) Or en : : ; : 5 : 
S Sagittae | 
SU Cygni 070 +0.25 —0.08 -—0.04 0.00 40.04 @ 
OFO .—0.07 =0.31 0/10 0:06 “041 047) Ort “Os 0s, oe 0.00 —0.01 —@ 
0.1 +0.03 —0.21 .—0.06 —0.02 +0:08 +0.10| 0.2 +0.43 +0.08 +0.04 0,00 se=vouue 
0.2 +£0.08 —0.13 =0.03 =0-01 4-08 10.03 141053" cove) y-LOvie) 70.07 0.00 -—0.08 —0. 
0:3: 402 2007 0.00 —0.01 +0.01 —0.04| 0.4 +0.77 40.31. +0.12 —0\01 70s een 
0.414016. 0.02 4003 0.00 —0.02 —0.07 | 0:5 +0.95 40.39 40.15 —0.01 =09aume 
0.5-. 40.18 -+-0.02 +40.04 0.00 —0.03 ~007| 0.6 +1.09 +0.45 -+0.18 —0,01) one 
0.6 +0.18 +0.04 +0.05 —0.01 —0.03  -—0.07 | 0.7 +1.05 +0.41' +017 —0.01) 0 ee 
0.7. +0.14 +0.01 40.04 —0.02 -—0.01 —0.05| 0.8 +0.68 40.27. +0.10 | 0.00 =m 
0.8 0.00 ~—0.12° —0:01 90:02. 40.04 440.07 | (0;9le 0.84 | oe 0.00 —0.02 —Om 
0:9'< <OAS A034 2042 001 wsOs124 GEO 
X Cygni 
T Vulpeculae 070 +0.91 +0.30 40.18 -0.01 -0.16 04 
070 —0.01 -0.27 -0.11° -0.01 40.13 +40.19| 0-1 aa a co mete 
0:1 —oi0g - 0:19 =0/08 "= 0.00;-.-0.08! 0,04 : ; ‘ : : ; r* 
2 : 0:3 +1.27  +0.64 +0.29° —0:02) )=0:3ene eam 
0.2 +0.19 0.10 —0.03 —0.01 ‘40:04 -+0.07 
0:3 40:29 0:01 40:01 —0.02* 0.00. 0.00 | 0-4 +1.66 +0.78 +0.34  —0.02) )==0-siiiamm 
0.4 40:38 40:07 40:04 —0.01 —0.03 —0.06| 0-3 71:82 +0.88 -+0.38 | —0.02 0-0 ges 
0:5 41044 40.13 40.06. 0.01 ° 0.04 0.10. 0-6 1.81 0.93" “10:41 ) 0:02 —seamas =o 
; : i é ; 0.7 +1.81 +0.97 +0.40 —0.01 =010nu—emm 
0:6 +048 50.96, “006° 0:01 =-0:02 0109 y 
‘a ee 0.8 +1.75 +0.97 -+0.36 .--0:01\ Gre nee 
O27. 0-460 3 016) 450060 0001 0.05 —0.06 
0:8 40.35 40:06 10:02 “0.01 0.01 20.01 | 0-85 41-15 +0.70 +0.31 +0.027)=30 ssnni sam 
e ri . >: P . x 
0:9 40.09 017 0.05 0.01 40.07 40,13 | 0-9, +1-20 +0.54 0.30) --0.01 1) 70ss0sii osm 
: : : 0.94 41.23 +0.75 -40.29 -++-0:03  ~=0)32ni=ame 
FF Aquilae Y Ophiuchi 
0FO 40.39 0.00 +0.02 -—0.01 -—0.01 —0.12 |: OFO +1.36 +0.76 -+0.45- +0103) =077 eee 
Olan ROA 7a 4-004 7) 40.05 0:00 =0.05. —0.17} 0.1 +154 +0.84 +0.43 +2002) Oulu 
0.2 +0.54 +0.11° +0.08 0.00 —0:07  =—0.20 | 0.2 /41:70 +0.90 40.47 +-0.02)00= gman 
0.3 +0.59 +0.14 +0.09 0.00 —0.08 \—0.23| 0.3 +1.84 40.95 +0.51° 40.02) =0)3euneomm 
OA 0:62) £0.16 0111 0.00 —0.12 —0.25| 0.4 +12.95 +1.02 +0.55 © 0:02) == 095s ume 
0.5 40563). 40.19 <40-12 0.00 —0.13 -—0.26| 0.5 +2.00 +4.11 +40.56 +0,02 (= Osseuemn 
O16) =H0.61 4-06.10) 220. 13 0.00 —0.13 —0.27 | G16 +2:00 +1.15 +0.57 +0.01) = 0)s0uee 
0.7 .4+0.58 40.17 +0.10 0.00 —0.10 —0.24] 0,7 +1.95 +1.08 40.56 10.01 = =00ecuueee 
0:8 £0.51. -4-0514- £0.07 0.00 —0:.08 —0.19 | 0.8 +1.83 +0.99 -+0.54' +002) 2075 queen 
0.9 +0.45 +0.05 +0.05 =0:01 -—0.05 -—0,14} 0.9 +1.63 -++0.89 £0.50 -10/0300=-G0sgee 


hydrogen absorption during the maximum phase. The 
deviation of V is positive and is a function of the spectral 
type at minimum light as given by Code (1947), as well 
as of the range of the spectral type from maximum to 
minimum. It is probably due to the presence of nu- 
merous metallic absorption lines in this spectral region, 


which increase in strength at minimum light. The regio 
of the I color is near the maximum of the negativ 
hydrogen-ion absorption, and the positive deviation 0 
the I point may be due chiefly to changes in that ab 


sorption. 


If we derive from the above-mentioned plots th 


CEE AE 912 HiO a ONE Rey 


04 0.8 


12 Log P 


1G. 3. Slope of the range-reciprocal wavelength relation as a 
ction of the logarithm of the period. Spots: type A and B 
heids; circles: type C cepheids; cross: SU Cassiopeiae. 


pes of the straight lines, using only the three colors B, 
and R, and plot these values, which we shall call 
s and which are given in Table VII, as a function of 
> logarithm of the period, we find that the various 
ints define two straight lines, well separated from 
+h other, as can be seen in Fig. 3. This diagram also 
ludes points for 6 Cep (type A), 7 Aql (type B), 
JMi (type C) observed by Stebbins (1953 and 1946), 
1 for SZ Tau and ¢ Gem (both type C), taken from 
ent work by the writer (Svolopoulos 1959), The 
per line has the equation k=1.00-logP—0.175, and 
ates to A and B-type cepheids, while the lower has 
equation k=0.275-logP+0.075, and relates to type 
sepheids. The 15 stars plotted in Fig. 3 are scarcely 
squate to establish such a separation for all cepheids; 
yertheless, the scatter in such a plot should be less 
in that in the plots attempting to separate the various 
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TABLE VI. Phase retardation of light maximum between U and I. 


SU Cas 0703 Y Ser 0.00 
DT Cyg 0.05 U Aql 0.06 
SU Cyg 0.01 S Sge 0.03 

T Vul 0.04 X Cyg 0.02 
FF Aql 0.02 Y Oph 0.06 


types of cepheids by use of spectral regions that are 
more affected by hydrogen or metallic-line absorption. 
For example, the plot given by Eggen (1951) relating 
photographic range to period will be more affected by 
such absorption; this effect is shown by the deviations 
of our V points from a linear plot. The physical inter- 
pretation of the result shown in Fig. 3 is not clear to: 
the writer. 

We have attempted to compute the values of k 
directly from our color data in a different manner. For 
a blackbody, the Wien approximation gives for the light 
range, Ay, of a variable star a 1956): 


1 max 
JA) log Sees) a ei a a 
min er iia 


where 7max and 7min are the maximum and minimum 
radii and Tmax and Tmin are the corresponding surface 
temperatures. We then have 


Ges) 
Reomp = — € rate IC 
: I races 1 satin 


The color temperatures of the various cepheids at both 
maximum and minimum light were computed directly 
from the observed six-color data in the following manner. 
The B-R color differences, which refer to a range of 
spectrum not much affected by absorption lines and 
bands, were first corrected for color excess, and then 
the corresponding blackbody temperatures were derived 
from the table given by Stebbins and Whitford (1945). 


TABLE VII. Temperature data for fifteen cepheids. 


Cepheid 

Star type CE (p_r) (B—R)max (B—R)min Trax Sets Robs Room 
SU Cas ? +0"37 —0.30 —0.12 6650 5900 0.21 0.30 
DT Cyg C +0.13 —0.35 —0.20 6900 6200 0.20 0.26 
SZ Tau G +0.48 —0.28 —0.07: 6550 5900: 0.25 0.26 
SU Cyg A +0.30 —0.49 —0.20 7650 6200 0.50 0.48 
a UMi (e +0.10 —0.26 —0.20 6750 6200 0.10 0.21 
T Vul A +0.23 —0.46 —0.11 7500 5850 0.49 0.59 
FF Aql Cc +0.40 —0.36 —0.16 6950 6050 0.27 0.33 
6 Cep A +0.22 —0.51 —0.02 7800 5600 N72 0.79 
Y Sgr A +0.34 —0.30 +0.12 6650 $150 0.65 0.69 
U Aql B +0.51 —0.17 +0.16 6100 5000 0.59 0.56 
n Aql B +0.30 —0.37 +0.06 7000 5300 0.64 0.72 
S Sge B +0.20 —0.23 +0.15 6300 5050 0.71 0.61 
¢ Gem Cc +0.17 —0.16 +0.09 6000 5200 0.33 0.49 
X Cyg B +0.41 —0.17 +0.39 6100 4550 1 il 0.87 
Y Oph Cc +0.85 +0.03 +0.29 5600 4750 0.43 0.49 
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TaBLeE VIII. Radius variations for four cepheids. 


dr/dt r—R dr/dt r—R 
Phase (km/sec) (105 km) Phase (km/sec) (105 km) 
DT Cygni T Vulpeculae 
0?0 +9.6 —1.7 0F0 +25.3 —8.9 
Ord +11.2 +0.6 0.1 +20.2 +0.0 
0.2 +5.4 42.3 0.2 +12.7 +6.3 
0.3 +1.0 +3.0 0.3 +5.5 +9.9 
0.4 —3.0 +2.8 0.4 —1.8 +10.6 
0.5 —6.1 +1.8 0.5 —8.9 +8.5 
0.6 —8.2 +0.2 0.6 —15.5 +3.8 
0.7 —10.7 —1.7 0.7 —20.6 —3.3 
0.8 —4.9 —3.5 0.8 —19.8 —11.2 
0.9 +5.1 —3.5 0.9 +3.8 —15.6 
SU Cygni Y Ophiuchi 
0r0 +24.7 —4.5 0r0 +5.1 —23.0 
0.1 +20.5 +3,.2 0.1 +10.3 —10.0 
0.2 +12.0 +8.6 0.2 +12.0 +6.0 
0.3 +3.5 +11.2 0.3 +7.5 +21.0 
0.4 —4.9 +11.0 0.4 +0.3 +27.0 
0.5 —16.2 +7.5 0.5 —5.8 +23.0 
0.6 —27.5 +0.2 0.6 —10.7 +10.0 
0.7 —27.5 —9.6 0.7 —11.1 —7.0 
0.8 —4.9 —15.6 0.8 —6.9 —20.0 
0.9 +23.3 —12.5 0.9 —1.3 —26.3 


The color excesses in B-R were derived from the P-V 
color excesses given in a recent paper by Kron and 
Svolopoulos (1959), in the same way as those in P-V 
were derived from those in V-B by Kron (1958b). The 
color excesses, B-R values, derived temperatures, and 
computed values of & are given in Table VII. The com- 
puted values of & are, generally, somewhat larger than 
the observed values, especially for the type C cepheids, 
but in view of the uncertainties in the computation the 
agreement seems fairly satisfactory. The temperatures 
and computed & values are, of course, quite dependent 
on the interstellar absorption values adopted. 

The temperatures of the type C cepheids'appear to 
be somewhat lower at maximum light and somewhat 
higher at minimum light than those of the type A and 
B cepheids of similar period. The temperature ranges 
are much larger for the A- and B-type cepheids than 
for those of C-type, as anticipated: AT4—3=1550°, 
ATc=750°; the mean temperatures of the three types 
of cepheids at the same period do not appear to differ 
appreciably, within the errors of the present deter- 
minations. 

Combining the relations giving the light range as a 
function of inverse wavelength and the straight line of 
Fig. 3 for cepheids of types A and B ,we obtain 


Am= (logP—0.175) -w4—1.05 logP-++0,42, 


which gives the light range as a function of the period 
and the wavelength. A check of this formula with the 
observational data for the colors B, G, and R was quite 
satisfactory for all of the cepheids discussed here. Also 
for U—for those stars unaffected by absorption in that 
color—the formula was also satisfactory. Since no stars 
of period around 9 days were measured, it was not 
possible to test a possible discontinuity of the relation 
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there, as others have suspected (Ledoux and Walra 
1958). As only a few type C variables are includec 
our list, such a relation was not derived for them. 
For four of the cepheids that have reliable rac 
velocity curves, an attempt was made to determine 
mean radius by applying the Wesselink (1946) mett 
in the same way as Stebbins (1953) had appliec 
previously. The rate of variation of the cepheid’s 
was computed from the radial-velocity curve assum 
a reduction factor of 24/17, which is appropriate 
intermediate limb darkening. The radial-velocity cut 
used were corrected in phase, utilizing the element 
the present photometric observations. Numerical i 
gration then gave the radius displacements, wh 
when combined with the difference in magnit 
between two phases of equal color, gave the m 
radius. The variation of the radius, dr/dt, and 
radius displacements, y—R, are given in Table VIII 
the four cepheids. The displacement curves are plot 


o’e of0 Of2 04 O% ofa 1%0 


Fic. 4. Radius displacement curves for four cepheids. 
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TasBLeE LX. Computation of cepheid radii. 


_DT Cygni 

Phase 1 0F05 OP15 0?40 OF45 
Phase 2 0P895 OF85 0°63 0?60 
Am 0060 07068 07072 07050 
n 1.0280 1.0318 1.0337 1.0239 
Ari: 10-5 km —0.5 +1.6 +2.7 +2.3 
Ar2:10~-> km —3.6 —3.8 —0.4 +0.2 
R-1077 km iba fol Dat 0.9 
Weight 1 3 2 Ry 

SU Cygni 

Phase 1 0F00 0?20 0?30 

Phase 2 0F87 0P805 OP77 
Am 0™075 0.116 0.092 
n 1.035 1.055 1.0435 
Ar;:10~> km —4.5 +8.5 +11.2 
Arz:10-5 km —14.7 —15.6 —14.8 
R-1077 km 3.06 4.5 6.13 
Weight 1 4 2 

T Vulpeculae 
Phase 1 Osalbys 0F25 0F35 0?40 OF45 
Phase 2 0P89 OF85 0P815 0?795 OP775 
Am 07176 07144 0™122 0™112 07094. 
nN 1.0845 1.0686 1.058 1.053 1,044 
Ar,:107* km =++3.5 +8.4 +10.5 +10.6 +9.8 
Ar2:10-5 km —15.6 —14.6 —12.3 —10.8 —9,2 
R-1077 km 2.4 340: 4.1 4.2 4.4 
Weight D 3 3 3 2 

Y Ophiuchi 
Phase 1 OP15 0°20 0F30 0F40 
Phase 2 0P9o1 0?90 OF85 OP77 
Am 0=252 0=186 0™158 07130 
n 1.123 17090 1.075 1.062 
Ar,:10~§ km —1.0 +6.0 +21.0 +27.0 
Ar2-10~5 km —26.5 —26.4 —24.0 —16.3 
R:1077 km 2.34 3.81 6.40 7.10 
Weight 1 1 1 il 


ig. 4. Table IX gives the pairs of phases of equal 
-in the ascending and the descending branches 
were matched, and also details of the computation 
e radius by the formula 


Ari—nAre 
R=-———., 
n—1 

e Ar;, and Arp are the radius displacements at the 
matched phases 1 and 2, and 7 is computed from 
elation logn=0.2 Am, where Am is the mean dif- 
ce in magnitude in the four colors V, B, G, and R. 
the matched phases 1 and 2 are not at the extreme 
ts of maximum and minimum light, but between 
, the deviation in Am of V, because of the metallic- 
ibsorption, is smaller than the observational error, 
ere is no reason not to use the V color in this radius 
utation. For the same reason the colors U and I 

occasionally used.) Table X gives the mean, 
mum, and minimum radii of the stars, and the 
ional variation in the radius. The probable errors 
1e weighted means of the individual values of R 
so given in the table. For Y Ophiuchi the results 
n in Table X agree very well with the computation 
bt (1954), who used Eggen’s two-color photometric 
rvations, 


For SU Cassiopeiae the application of the Wesselink 
method led to unsatisfactory results. SU Cassiopeiae is 
bluer in a part of the descending branch near minimum 
light than it is in the ascending branch at points where 
the brightness, as observed in the four colors used, is 
the same; the computation thus gives a negative radius 
near minimum light and the mean radius is indeter- 
minate by this method. It would seem that this star 
must deviate considerably from a blackbody, as in cases 
previously discussed by Arp (1957). 

From the radii and the radius-displacement curves, 
the light variation produced by the changing surface 
area of each of the four stars was computed. Subtracting 
the surface-area effect from the mean observed light 


TABLE X. Summary of radius data. 


Rmax Rosin 
(107 (107 Reference for radial- 
Star (R/Ro)mean km) km) dR/R velocity curve 
DT Cyg 2243 £1.53 1.46 0.05 sari te and Herbig 
1952) 
SUCyg 7048 5.0 4.7 0.06 Becker and Strohmeier 
(1942) 
TVul 554+3 3.9 3.6 0.07 Liist-Kulka (1954) 
YOph 70+10 5.2 4.6 0.11 Abt (1954) 
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Oo +0.2 
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Fic. 5. Surface brightness in V as a function of the color V. 


curves, we obtain the variation of the surface brightness 
with phase. A plot of the surface brightness in the color 
V against the V color for the various phases gave a 
straight line for each star included in Table X (Fig. 5). 
This result is consistent with the assumption of the 
Wesselink method that the surface brightness is a 


SVOLOPOULOS { a 


single-valued function of the color. The effect 
metallic-line absorption on the V magnitude is 
small to appear in Fig. 5. 
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